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Abstract 

The  United  States  Air  Force  Research  Lab,  Munitions  Directorate,  Flight 
Vehicles,  Integration  Branch  (AFRL/MNAV),  has  developed  a  flexible  wing  micro  air 
vehicle  (MAV)  to  be  used  with  special  tactics  groups.  In  keeping  with  the  requirement 
that  this  MAV  be  compact,  previous  research  developed  a  rotatable  tail  mechanism  which 
resulted  in  tail  movement  similar  to  that  of  a  bird’s  tail.  In  this  study  the  design  of  this 
tail  was  modified  to  produce  a  more  storable  vehicle.  The  redesign  also  allowed  the  tail 
to  deflect  upward  to  large  angles,  enabling  the  tail  to  be  used  as  a  spoiler.  The 
aerodynamic  affects  of  adding  a  vertical  stabilizer  mechanism  to  improve  the  stability  of 
the  vehicle  and  rotatable  tail  combination  was  also  quantified.  Data  from  these  tests 
confirmed  the  tail  is  a  plausible  method  to  reduce  lift  and  increase  drag,  consistent  with 
proper  spoiler  function.  A  wide  range  of  angles  were  used  to  demonstrate  that  forces  and 
moments  from  the  flat  tail  were  similar  to  those  of  traditional  rudder.  Directional 
stability  was  improved  by  the  stabilizer,  and  recommendations  for  further  improvements 
are  given. 
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CHARACTERIZATION  OF  A  ROTARY  FLAT  TAIL  AS  A  SPOILER  AND 
PARAMETRIC  ANALYSIS  OF  IMPROVING  DIRECTIONAL  STABILITY  IN  A 

PORTABLE  UAV 


I.  Introduction 


Background 

In  today’s  ever  changing  battlefields  enemies  are  using  more  unconventional 
means  of  attack  and  evasion.  This  can  be  seen  most  clearly  in  the  current  global  war  on 
terrorism.  Increasingly  unorthodox  approaches  must  be  developed  to  seek  out,  monitor, 
and  attack  these  enemies. 

The  United  States  Anned  Forces  currently  employs  the  Predator  and  Global  Hawk  for  its 
long  range  reconnaissance  work.  The  U.S.  also  is  developing  and  using  smaller  UAVs 
(Unmanned  Aerial  Vehicles)  to  fill  the  need  for  reconnaissance  and  battle  damage 
assessment  on  the  company  or  platoon  level  (Parga,  2004:  3).  The  current  UAVs  for  this 
type  of  work  are:  Desert  Hawk,  Shadow  and  Dragon  Eye.  Room  does  exist  for 
improvements  in  size  and  portability  of  these  UAVs. 

The  United  States  Air  Force  Research  Lab,  Munitions  Directorate,  Flight 
Vehicles,  Integration  Branch  (AFRL/MNAV)  has  produced  a  flexible  wing,  portable 
UAV  to  study  and  develop  a  better,  smaller  unit  UAV  for  the  war  on  terrorism.  This 
vehicle  was  studied  by  Captain  DeLuca  and  Lt.  Rivera  Parga  in  their  theses.  The  focus  of 
this  study  is  to  build  upon  the  “bird  like”,  flat,  rotating  tail  described  by  Lt.  Rivera  Parga. 
The  influence  of  avian  research  is  not  limited  to  only  Lt.  Parga ’s  work;  this  research  may 
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prove  to  be  beneficial  to  the  development  of  future  UAVs.  Birds  demonstrate 
considerable  tail  muscle  control,  which  gives  them  considerably  more  degrees  of  freedom 
than  the  two  degrees  used  in  this  and  previous  studies. 

Today,  avian  flight  is  used  to  improve  the  understanding  of  flight.  Avian  tails  are 
being  researched  to  find  different  manners  to  control  flight  and  improve  upon  current  tail 
designs.  Studies  conducted  by  R.G.  Hoey  and  then  A.  L.  Thomas  suggest  several 
possibilities  of  how  birds  use  their  tails  and  how  this  may  develop  into  use  by  humans. 
Their  research  has  uncovered  several  manners  of  which  birds  use  their  tails. 

However,  “The  study  of  bird  tails  is  very  polarized.  One  camp  of  theorists 
emphasizes  the  importance  of  ordinary  physics,  and  the  other  camp  points  to  the 
extraordinary  tails  of  male  peacocks  and  barn  swallows  as  examples  of  sexy  fashions 
overpowering  sensible  aerodynamics”  (Rayner;  1988). 

Examples  that  the  researchers  who  favor  the  use  of  bird’s  tails  for  control 
offer  strike  a  resemblance  to  the  use  of  tails  on  aircraft.  For  instance,  one  author 
(Horton-Smith  in  a  1938  text,  The  Flight  of  Birds,  page  38)  states:  “A  bird,  like 
an  airplane,  uses  rotation  of  wing  and  tail.  There  is  no  vertical  fin  in  the  bird’s  tail 
so  it  has  to  rely  on  banking.  It  is  possible  that  a  long  tail,  when  bent  to  one  side, 
may  function  as  a  rudder”.  Additionally,  John  H.  Storer  comments:  “The  tail  of  a 
bird,  indeed,  has  many  uses.  It  can  steer  in  any  direction,  act  as  a  brake,  form  a 
slot  behind  the  wings,  or  become  a  part  of  the  bird’s  lifting  surface, 
supplementing  the  wings.  The  swallow-tailed  kite  twists  its  tail  to  steer.  It  may 
turn  its  tail  so  that  either  the  upper  or  the  lower  surfaces  will  strike  the  air  stream 
in  steering.  The  sides  of  the  tail  may  be  controlled  separately.” 
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Opposing  theorists  offer:  ““The  tail  of  birds  virtually  has  no  stabilizing 
effect.  It  is,  hence  not  a  stabilizing  instrument.  Only  to  a  limited  extent  is  it  used 
as  a  steering  device.  Mostly  it  is  used  at  low  speed  as  a  landing  flap”  (Nickel  & 
Wohlfahrt,  1994:  25).  These  same  theorists  state  that  a  bird  achieves  stability  by 
varying  lift  generated  by  each  wing,  through  flexure,  changing  angle  of  attack 
independently,  and  sweep.  Warrick  et  al.  sum  the  debate  of  avian  tail  use  best 
with  the  following  statement:  “Beyond  its  theoretical  capabilities,  the  precise  use 
of  the  tail  in  flying  birds  has  not  been  thoroughly  documented.” 

However,  Warrick  et  al.  then  go  on  to  present  a  histogram,  Figure  1,  of 
occurrences  of  tail  elevation,  twist,  depression,  combined  twist  and  depression, 
and  no  tail  manipulation  for  a  swallow  during  a  series  of  videotaped  prey-capture 
maneuvers.  Though  the  data  is  limited  to  a  few  hundred  maneuvering 
observations,  the  data  suggests  that  tail  elevation  (pitch  up)  is  used  only  to 
increase  drag  in  order  to  stall  the  flight  and  initiate  a  dive.  Furthermore,  the  most 
common  tail  orientation  during  a  level  turn  was  a  combined  tail  depression  and 
twist  while  the  most  common  tail  orientation  during  a  climb  was  depression  (pitch 
down)  only.  By  and  large,  this  set  of  observations  is  consistent  with  comments 
presented  above  which  favored  tail  use  as  a  flight  control  mechanism. 


3 


Prt)  cijmireM.rirtivei 


Figure  1.  Uses  of  Bird  Tails  (Warrick  et  al,  2002) 

In  the  thesis  of  Jose  Parga,  micro  air  vehicle  tails  were  modeled  after  studying  the 
tails  of  various  birds.  Since  modeling  a  complete  bird’s  tail  with  its  many  movements 
would  be  almost  impossible,  a  design  that  resembled  a  bird’s  tail  in  combined  rotation 
was  deemed  more  practical.  The  result  was  a  bird-like  tail  that  could  be  used  to  affect 
pitch  as  well  as  roll.  However,  problems  with  stability  were  encountered  with  the  MAV 
and  vertical-stabilizer-less  configuration  that  birds  do  not  seem  to  display.  In  the  thesis 
by  Parga  it  was  summarized  that  the  micro  air  vehicle  in  question  would  be  able  to 
achieve  controllable  flight  if  modifications  to  the  tail  were  made. 


Motivation 

In  keeping  with  an  inherent  goal  of  micro  UAV  design,  building  a  smaller 
vehicle,  Lt.  Parga  removed  the  V-tail  configuration  on  the  AFRL/MNAV  micro  air 
vehicle  (MAV).  It  was  then  replaced  with  a  “rotary”  flat  tail  that  was  developed  after 
studying  the  tails  of  various  birds.  This  rotary  tail  enabled  the  tail  to  move  as  if  it  was  an 
elevator  and  at  the  same  time  revolve  about  the  x  axis  of  the  aircraft  body  (Parga,  2004: 
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10).  Having  characterized  this  tail  as  a  plausible  means  to  control  flight,  Parga’s  thesis 
then  summarized  control  effects  by  measuring  changes  in  the  force  and  moment 
coefficients.  Motivation  for  the  current  thesis  was  found  in  searching  for  ways  the 
MAV  could  become  more  recoverable.  One  way  to  accomplish  this  is  to  reduce  the 
horizontal  velocity  of  the  vehicle  prior  to  landing  by  using  the  tail  as  an  airbrake. 

To  accomplish  this  task  the  tail  was  envisioned  to  deflect,  as  an  elevator,  to  -90°. 
The  tail  used  by  Parga  could  not  accomplish  this  deflection  without  reducing  the  range  of 
positive  deflection.  A  different  tail  using  the  same  volume  coefficient  as  Lt.  Parga’s  tail 
designs  was  constructed.  The  new  tail,  however,  was  completely  flat  and  attached  to  the 
MAV  in  a  slightly  different  way.  The  configuration  change  also  allowed  the  tail  to  fold 
completely  parallel  with  the  top  of  the  MAV,  thus  improving  its  compactness  for  storage 
or  transport,  witnessed  in  Figure  2. 


Figure  2.  Fully  Foldable  MAV  Tail 

Due  to  Cnp  being  negative,  as  evidenced  in  Lt.  Parga’s  thesis,  a  stabilizing 
mechanism  was  developed.  This  mechanism  consisted  of  two  ventral  fins  placed  on  the 
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lower  aft  portion  of  the  MAV  fuselage.  This  arrangement  gave  the  added  benefit  of 
providing  a  hard  point  to  attach  rear  wheels  for  landing  to  make  the  vehicle  more 
recoverable. 

Research  Objective 

a.  The  primary  focus  of  this  study  is  to  use  the  AFIT  open-circuit,  low-speed 
wind  tunnel  to  characterize  the  behavior  of  the  improved  flat  tail  and  stability 
modifications.  To  achieve  this  end  state  the  following  data  points  were  gathered: 

b.  Calculate  and  compare  the  lift,  drag,  and  side  force  coefficients,  Cl,  Cd,  and  Cs 
for  the  Parga  tail  vs.  the  flat  tail  in  a  high  deflection  condition. 

c.  Calculate  and  compare  the  force  coefficients,  moment  coefficients  and  stability 
derivatives  for  larger  values  of  elevation  and  rotation  angles  of  each  tail. 

d.  Calculate  and  compare  force  and  moment  coefficients  of  the  two  tail 
configurations  in  various  deflection  and  rotation  angles. 

e.  Add  vertical  fins  and  recalculate  the  stability  derivatives  of  the  two  tails  for  this 
configuration. 
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II.  Literature  Review 


Introduction 

Extensive  research  into  unmanned  aerial  vehicles  (UAVs)  during  the  past  few 
years  has  led  to  advancements  in  the  miniaturization  of  cameras,  data  transmitting  links, 
navigation  systems,  thrust  generator,  lifting  surfaces,  and  control  surfaces  (Gad-el-Hak, 
2001:  419).  Miniaturization  of  lifting  surfaces  has  brought  about  its  own  set  of 
challenges. 

With  chord  lengths  of  micro  air  vehicle  wings  becoming  smaller,  Reynolds 
numbers  also  become  small.  Most  MAVs  fall  with  in  the  Reynolds  number  range  104  - 
106  (Carmichael,  1981).  When  Reynolds  number  falls  in  the  range  of  104  -  103' 
challenges  are  presented  because  the  lifting  surfaces  experience  greater  susceptibility  to 
flow  separation  and  lower  efficiency  (Gad-el-Hak,  2001:  419).  The  MAV  in  this  study 
has  Reynolds  numbers  ranging  from  0.9  X  105  to  1.3  X  105  for  the  velocities  tested. 

Low  Reynolds  Number  Design 

“Within  104  <  Re  <  106  a  complex  flow  phenomena  takes  place  on  the  upper 
surface  of  the  airfoil.  Separation  of  the  flow  will  take  place  before  transition  to  turbulent 
flow.  Flow  in  this  range  of  Reynolds  numbers  is  laminar  and  can  support  only  a  very 
small  adverse  pressure  gradient  without  separation.  With  any  surface  roughness  on  the 
wings  pressure  gradients  become  inherent.  Thus,  the  laminar  flow  separates  and  forms  a 
free-shear  layer  that  quickly  transforms  to  turbulent  flow.  Reattachment  of  the  flow  will 
then  occur  if  sufficient  energy  is  supplied  to  the  wall  region  by  entrainment  of  the  high¬ 
speed  fluid.  A  separation  bubble  will  form”  (Gad-el-Hak,  2001 :  419). 
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Figure  3.  Laminar  Separation  Bubble  (Gad-el-Hak,  2001:  421) 


A  separation  bubble,  shown  in  Figure  3,  will  serve  to  de-energize  the  flow  over 
the  wing  and  result  in  decreased  lifting  efficiency.  Aerodynamic  efficiency  of  the  wing 
is  already  low  within  the  range  of  Reynolds  numbers  due  to  viscous  forces  producing  a 
high  drag  in  the  shear  layer.  According  to  Nechyba  and  Ifj'u,  biological  aircraft  such  as 
bats,  insects,  and  birds  are  able  to  mitigate  this  problem  and  are  able  to  outperform  any 
manmade  miniature  flight  vehicle. 

The  optimal  situation  for  MAV  designers  would  be  to  create  a  scenario  where  the 
size  of  the  separation  bubble  decreases  as  the  incident  angle  of  attack  increases  through 
the  stall  angle,  astaii,  at  which  point  the  flow  enters  a  turbulent  transition  state  near  the 
trailing  edge  (Biber  et  al,  2004:7).  Using  an  aeroelastic  wing  with  the  ability  to  adapt  to 
atmospheric  disturbances;  researchers  with  The  University  of  Florida  and  NASA  Langley 
Research  Center  have  accomplished  such  a  task  (Waszak,  2001 :  1). 
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Flexible  Wing  MA  Vs 

The  University  of  Florida  research  group  is  considered  to  be  at  the  forefront  of 
MAV  design  with  root  chords  of  6”  or  smaller.  An  example  of  one  of  the  group’s 
vehicles  is  presented  in  Figure  4.  The  group  has  experienced  great  success,  winning  the 
International  Society  of  Structural  and  Multidisciplinary  Optimization  (ISSO)  Micro 
Aerial  Vehicle  Competition  for  three  years  in  a  row  (Waszak  and  Jenkins,  2001 : 1).  They 
accomplished  the  task  by  using  carbon  fiber  matrix  ribs  supporting  a  parachute-latex 
membrane  to  make  the  planar  surface  of  the  wing;  enabling  what  is  tenned  as  adaptive 
washout. 


Figure  4.  University  of  Florida’s  3.3"  Mean  Chord,  Flexible  Wing  MAV 


In  adaptive  washout  the  flexible  wings  may  compensate  for  changes  in  speed  and 
altitude  by  changing  angle  of  attack  along  the  span.  This  adaptive  washout  is  also  seen  in 
sails  of  sail  boats.  Here  the  sail  is  twisted  by  wind  velocity  to  extend  the  wind  range  of 
the  sail  and  produce  a  more  constant  thrust,  even  in  gusty  conditions  (Ifju  et.  al,  2002:  2) 
Local  flow  conditions  can  be  compensated  for  by  the  flexible  membrane  moving 
or  twisting  to  achieve  a  more  constant  lift  (Waszak  and  Jenkins,  2001 :2).  In  other  words 
the  severity  of  negative  pressure  gradients  is  reduced.  Wind  tunnel  tests  were  conducted 
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by  The  University  of  Florida  to  compare  differences  in  a  flexible  wing  and  rigid  wing 
design.  Figure  5  shows  the  improvement  in  lifting  characteristics  of  a  flexible  wing 
design.  Additionally  the  work  of  Capt.  Deluca,  Figure  6  and  Figure  7,  is  included  for 
comparison.  The  plots  show  that  the  flexible  wing  is  able  to  achieve  higher  L/D  and  lift 
coefficient.  Cl  plots  show  how  a  higher  angle  of  attack  can  be  achieved  with  the  flexible 
wing,  demonstrating  that  separation  of  flow  and  the  resulting  loss  of  lift  is  delayed. 


-m —  Rigid  (Graphite) 


—  Rigid  (Graphile) 


Figure  5.  L/D  and  CL  vs.  a  for  University  of  Florida  MAV  (Waszak  &  Jenkins,  2001:  4) 
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Figure  6.  L/D  Plot  for  Rigid  and  Flexible  Wing  MAY  (DeLuca,  2004:  54) 


Figure  7.  CL  Plot  for  Rigid  and  Flexible  Wing  MAY  (DeLuca,  2004:  48) 
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A  second  aspect  of  optimization  due  to  miniaturization  of  UAVs  is  the 
improvement  of  control  surfaces,  most  notably  the  tail.  The  tail  is  crucial  not  only  for 
control  but  for  stability.  While  some  aircraft  are  tailless,  most  use  a  conventional  T,  or 
cruciform  type  of  configuration  for  the  tail.  Small  glider  enthusiasts  have  been  using  a 
V-tail  configuration  due  to  it  compact  size  and  combination  of  control  surfaces.  The 
MAV  for  this  study  originally  had  a  V-tail  with  two  “ruddervators”.  This  configuration 
was  studied  by  Deluca. 

The  idea  of  a  rotatable  tail-plane  is  not  a  new  idea.  William  Nash  patented  a 
design  in  1992.  Patent  Number  5069143  is  assigned  to  a  tail  plane  that  can  rotate  from  a 
horizontal  stabilizer  to  a  vertical  stabilizer  position;  additionally,  it  may  rotate  about  an 
axis  parallel  to  the  y  body  axis,  as  a  regular  control  surface  would. 

In  the  same  year  as  the  patent,  Robert  G.  Hoey  conducted  research  into  avian 
flight  by  building  a  remote  control  raven  model.  The  objectives  of  his  research  were  to 
determine  if  soaring  birds  are  statically  stable  in  the  lateral  axis,  pitch  axis,  find  the 
source  of  this  stability,  and  to  investigate  the  method  birds  used  to  control  turns.  Two 
raven  models  were  built.  The  first  model  had  drag  flaps  on  the  wing  to  induce  yaw  and 
the  tail  was  deflected  as  an  elevator  only.  The  second  model  retained  the  drag  flaps  but 
the  tail  was  modified  to  allow  tail  rotation  and  elevator  deflection.  With  this  second 
model,  different  wing  plan  form  shapes  as  well  as  tip  feather  spread  and  its  affects  on 
stability  were  investigated.  In  R.G.  Hoey’s  paper,  “Research  on  the  Stability  and  Control 
of  Soaring  Birds”,  no  quantitative  analysis  of  the  tail  was  given.  Hoey  concluded,  “The 
handling  qualities  were  not  very  comfortable  for  a  human  pilot  but  they  are  probably 
completely  normal  to  a  raven.” 
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While  the  use  of  a  bird’s  tail  in  flight  is  highly  debated,  some  research  indicates 
that  tails  do  serve  a  purpose.  A  bird  bends  its  tail  to  one  side  using  it  as  a  rudder  (Horton- 
Smith,  1938;  38).  Proponents  of  birds  using  their  tails  for  control  present  evidence  that 
the  tail  may  be  used  as  a  lifting  surface  (Storer,  1948;  38).  With  any  distance  between 
the  tail  and  wings  of  a  bird,  the  tail  can  be  used  to  produce  a  moment,  thus  it  becomes 
similar  to  an  elevator  and/or  rudder. 
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III.  Methodology 


Introduction 

This  study  used  four  primary  pieces  of  equipment.  The  first  two  pieces  are  the 
MAV,  provided  by  AFRL/MNAV,  and  the  two  rapid  prototyped  tails.  The  AFIT  low 
speed  wind  tunnel,  the  third  piece  of  equipment,  provided  the  means  to  test  the  tails. 

Last,  the  AFIT-1  balance  was  used  to  measure  the  forces  of  the  MAV  in  the  wind  tunnel. 


MA  V  and  Tail  Description 

The  original  MAV,  shown  in  Figure  8,  is  constructed  of  a  carbon  fiber  matrix 
used  for  the  entire  fuselage  and  tail.  The  wings,  however,  use  the  carbon  fiber  for  the 
forward  %  of  the  chord.  From  this  %  chord  position,  carbon  fiber  ribs  are  then  used  to 
support  a  parachute  type  material.  (DeLuca,  2004: 16)  This  structure  is  rigid  in  upward 
bending,  but  is  flexible  in  downward  bending.  This  allows  the  wings  to  be  wrapped 
under  the  body,  increasing  compactness  for  storage.  Properties  of  the  original  MAV  and 
tail  are  given  in  Table  1. 


Table  L  Original  MAV  and  Tail  Dimensions 


Winq 

Root  Chord 

Span 

Area 

Aspect  Ratio 

Mean  Aerodynamic  Chord 
Leading  edge  thickness 
Parachute  Material  Thickness 

6" 

24" 

93.5  in2 
6.16 

4.2" 

0.025" 

0.005" 

Tail 

Mean  Chord 

Span 

Area 

Aspect  Ratio 
Thickness 

Horizontal  Tail  Volume  Coeff. 

2.35" 

6.3" 
14.8  in2 
2.7 

0.03" 

0.54 
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Figure  8.  Original  MAY  Selected  Dimensions  (Parga,  2004:  51) 


To  conduct  the  research  done  for  Parga’s  thesis  the  V-tail,  shown  in  Figure  8  was 
removed  and  a  rotatable  tail  mechanism,  Figure  9,  installed.  The  same  actuators  as  the 
V-tail  design  were  used,  although  they  were  repositioned.  Removing  the  old  tail  and 
attaching  the  new  tail  gave  an  overall  length  of  approximately  14.08  inches  vs.  18.2 
inches  of  the  original  design  wing  dimensions  remained  the  same.  When  folded,  the  full 
vehicle  length  was  less  than  10  inches. 
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Figure  9.  Foldable  Rotatable  Tail  MAV 

The  rotatable  mechanism  allows  the  tail  to  be  rotated  around  the  “x”  body  axis  for 
roll  control  as  well  as  to  be  deflected  for  pitch  control.  The  rotation  mechanism  was 
drawn  in  Solid  Works™  and  produced  on  a  Stratasys  3300  rapid  prototype  machine. 


Figure  10.  Rotation  Mechanism 


A  control  rod  attached  to  a  rack  runs  axially  through  the  piece  pictured  in  Figure 
10.  The  geared  rack  moves  fore  and  aft  along  the  x  body  axis,  depending  upon  the 
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control  command  given.  A  gear  positioned  on  the  rotation  mechanism  and  attached  to 
the  tail  is  meshed  with  the  rack.  As  the  rack,  seen  in  Figure  1 1,  is  actuated,  the  gear 
rotates  causing  a  deflection  in  the  tail.  The  control  hom  near  the  top  of  Figure  10  is 
attached  to  the  servo  that  controls  rotation.  The  control  horn’s  main  function  is  to  allow 
more  leverage  and  rotation  angle.  Actuating  the  servo  results  in  a  rotation  of  the  entire 
rotation  mechanism  and  tail.  Inside  the  fuselage  the  piece  in  Figure  10  is  installed  into  a 
yoke  that  constrains  all  movements  except  rotation  in  the  “x”  body  axis.  Figure  12  shows 
the  internal  configuration  and  control  actuation  attachments  of  the  MAY. 


Figure  11.  Rotatable  Tail  Attachment 
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Figure  12.  Internal  Layout  of  Rotatable  Tail 


The  MAV  tested  in  this  study  differs  slightly  from  the  one  in  Lt.  Parga’s  study.  A 
new  rotation  mechanism  was  printed  out  for  the  new  rapid  prototype  machine.  The  new 
machine  allows  for  tighter  tolerances  and  a  stronger  part.  A  second  benefit  of  the  new 
mechanism  is  the  allowance  for  increased  rotation  angles  over  the  one  tested  previously. 
This  was  desirable  as  this  study  sought  to  investigate  such  angles.  Additionally,  the 
servos  were  replaced  with  a  higher  torque  servo;  to  reduce  some  of  the  tail  movement  due 
to  wind  forces  during  tests. 

In  some  tests  conducted  in  this  study  vertical  stabilizers  were  attached  to  the 
model  to  improve  problems  with  directional  stability  as  uncovered  in  Lt.  Parga’s  work. 
The  stabilizer  was  placed  at  the  aft  end  of  the  aircraft,  just  forward  of  the  tail  mechanism. 
Two  screws  were  already  located  on  the  model  to  hold  the  tail  yoke  in  place.  The  screws 
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provided  an  excellent  attaching  point.  The  location  chosen  could  double  as  an  excellent 
position  for  landing  gear  supports.  To  examine  this,  two  holes  were  drilled  in  the 
stabilizers  and  landing  gear  axils  and  wheels  were  attached  for  later  tests. 


Figure  13.  Stabilizer 


The  stabilizers,  Figure  13,  were  designed  by  considering  the  second  tail 
configuration  in  Parga’s  thesis  which  had  vertical  stabilizers.  Results  from  Parga  showed 
that  the  vertical  stabilizers  improved  directional  stability  to  the  point  that  Cnp  was 
positive.  Stabilizer  area  was  increased  over  that  of  the  Parga  model  due  to  their 
placement  and  loss  of  moment  arm.  The  increase  in  this  area  was  weighted  against  the 
length  attainable  with  out  interfering  with  tail  operation. 


Table  2.  Stabilizer  Properties 


Length 

2.75” 

Height 

1.5” 

Width 

2.125” 

Area 

7.13  in2 

Weight 

0.040  kg 

The  tail  used  by  Lt.  Parga  and  labeled  as  “Tail  1”  in  his  tests  was  also  used  in  the 
tests  for  this  study.  It  will  be  referred  to  as  Tail  1  in  this  study.  The  use  of  this  tail 
provided  a  baseline  for  comparison  to  this  earlier  work.  Tail  1  was  designed  using  Solid 
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Works™,  printed  with  a  rapid  prototype  machine,  and  covered  with  a  polyester  fabric 
(Parga,  2004:  54).  Table  3  shows  the  properties  of  this  tail. 


Table  3.  Geometric  Properties  of  Tail  1 


Area 

9.42  in2 

Chord 

4.66" 

Span 

4.375" 

Thickness 

0.07" 

Aspect  Ratio 

2.02 

Taper  Ratio 

0.3755 

Tail  Volume  Coefficient 

0.2 

Parga  stated  that  the  total  estimated  time  for  construction  of  this  tail  to  be  9  hours. 
The  second  tail  used  in  this  study  consumed  a  similar  amount  of  time.  This  tail  will 
henceforth  be  referred  to  as  Tail  3  since  the  name  Tail  2  was  used  in  Parga’s  thesis  for  the 
tail  with  stabilizers.  Tail  2  is  similar  to  Tail  1  except  it  has  been  flattened,  whereas  Tail  1 
was  curved  laterally  as  evidenced  in  Figure  14  and  Figure  15. 


Figure  14.  Rear  Profile  of  Tail  1 


Figure  15.  Rear  Profile  of  Tail  3 


Tail  3  was  designed  using  Tail  1  as  a  guide.  Their  internal  structures  are  similar. 
Besides  being  flat,  Tail  3  differs  from  Tail  1  in  that  its  attachment  points  are  offset  from 
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the  plane  of  the  tail,  Figure  16.  This  offset  allows  Tail  3  to  be  fully  folded,  as  it  offsets 
the  tail. 


Figure  16.  Attachment  of  Tail  1  (top)  and  Tail  3  (bottom) 


Due  to  the  offset  in  attachment,  Tail  3  is  now  closer  to  the  vertical  position  of  the 
CG.  This  means  less  vertical  variation  in  CG  position  as  the  tail  is  deflected.  Also,  when 
rotated,  the  tail  now  rotates  around  the  x  body  axis;  before,  Tail  1  was  positioned  higher 
on  the  model  resulting  in  a  slight  lateral  CG  shift  when  it  was  rotated. 

Wind  Tunnel  Description 

The  AFIT  low  speed  wind  tunnel  is  located  in  building  644,  room  LI 54.  It  was 
constructed  by  New  York  Blower  Company,  who  also  manufactured  the  fan.  The  tunnel 
is  an  open  circuit  design  that  draws  and  ejects  air  from  within  the  room.  The  test  section 
is  closed  and  measures  approximately  3’  x  3’.  This  tunnel  has  a  design  test  speed  of  150 
mph  and  has  been  tested  to  a  speed  of  148  mph.  The  low  speed  wind  tunnel  can  be 
broken  down  into  6  components:  the  inlet,  converging  section,  test  section,  diverging 
section,  exhaust,  and  fan.  Figure  17  depicts  each  of  these  sections. 
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The  inlet  is  sized  to  dimensions  of  1 1 1”  high  x  122”  wide  x  70”  deep.  In  the  front 
of  the  inlet,  where  room  air  is  first  introduced  into  the  tunnel,  a  quarter  inch  aluminum 
honeycomb  structure  acts  as  a  flow  straightener.  This  is  followed  by  four  anti-turbulence 
screens.  The  idea  of  the  large  inlet  section  is  to  reduce  velocity  gradients  of  the  air 
entering  the  tunnel.  The  average  turbulence  intensity  of  the  tunnel  has  been  measured  to 
be  2.25%  (DeLuca,  2004:  85). 

After  the  inlet  section  the  air  is  channeled  and  sped  up  in  the  converging  section. 
Here  the  tunnel  converges  from  a  height  of  11”  to  a  height  of  3 1 .5”;  giving  the  section  a 
contraction  ratio  of  9.5:1.  The  converging  section  leads  into  the  test  section. 

The  test  section  measures  31”  high  x  44”  wide  x  72”  in  length.  The  top  and  sides 
are  constructed  of  Plexiglas.  The  sides  function  as  moveable  doors  for  accessing  the  test 
section.  The  top  consist  of  a  removable  Plexiglas  panel  with  openings  for  a  traversing 
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hot-wire  anemometry  system.  Within  the  test  section  is  the  sting,  used  for  mounting  the 
balance.  In  turn  the  sting  is  mounted  on  a  support  system  attached  to  a  turn  table. 

The  sting  support  system  is  able  to  move  20  degrees  up  or  down  from  its  neutral 
zero  angle  of  attack  position.  The  turn  table  is  able  to  rotate  from  sideslip  angles  of  +  15 
degrees  to  -15  degrees.  These  movements  are  controlled  from  the  control  room,  located 
in  room  LI  54  using  Lab  View  Virtual  Instrument®  software.  The  user  interface  of  this 
software  is  presented  in  Figure  18. 


Figure  18.  Wind  Tunnel  Control  Interface  (Parga,  2004:  75) 

b  24" 

The  wing  span  to  tunnel  width  ratio  is:  —  = - =  .0545  «  0.55  (DeLuca,  2004: 

w  44" 

25).  The  general  rule  is  that  wing  span  to  tunnel  width  should  be  —  <  0.8  (Barlow  et  ah, 

w 

1999:28).  Having  a  ratio  less  than  0.8  will  reduce  large  interferences  between  the  sides 
of  the  model  and  shed  vortices  from  the  wing  (Barlow  et  al,  1999:  381). 
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After  moving  through  the  test  section,  air  flows  to  the  diverging  section.  In  this 
section  air  is  diffused  to  reduce  its  speed  without  flow  separation  occurring  (Barlow  et  al, 
1999:  80),  The  length  of  the  diffuser  section  in  the  AFIT  low  speed  wind  tunnel  is  3 12”. 
The  diffuser  connects  to  the  fan. 

The  fan  is  an  ACFL/PLR  Class  IV  type  fan  manufactured  by  New  York  Blower 
Company.  It  is  powered  by  a  Premium  Efficiency  electrical  motor  manufactured  by 
Toshiba.  This  in  turn  is  controlled  by  an  adjustable  frequency  controller  manufactured  by 
Siemens.  Table  4  shows  motor  characteristics.  Upon  exiting  the  fan,  air  is  exhausted 
upward  into  room  LI 54. 


Table  4.  Wind  Tunnel  Motor  Properties  (Parga,  2004:  68) 


Type 

EQP  III,  3face  induction 

Power 

200  BHP 

Poles 

4 

Max.  operating  speed 

1785  RPM 

Voltage 

230/460  Volts 

Frequency 

60  Hz  Max 

Amperage 

444/222amp. 

Balance  Description 

AFIT-1,  a  six-component  internal  balance  was  used  in  tests  conducted  for  this 
thesis.  This  balance  was  manufactured  by  Modern  Machine  &  Tool  Company,  Inc. 


24 


Figure  19.  Position  of  Strain  Gages  in  AFIT-1  Balance  (Modern  Machine  &  Tool  Co.,  2004) 


For  the  tests  conducted,  the  balance  was  placed  inside  a  mounting  block  and 
attached  to  the  model.  The  mounting  block  is  geometrically  similar  to  a  camera  pod  that 
would  be  used  on  the  MAV,  and  the  intrusive  effects  are  not  considered  here. 

The  balance  is  composed  of  strain  gauges.  Each  strain  gauge  is  composed  of  a 
thin  wire  filament  wound  in  a  serpentine  pattern.  The  gauges  are  placed  in  a  Wheatstone 
bridge,  and  a  voltage  is  supplied  across  the  bridge  continuously  (Parga,  2004:71). 
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Figure  20.  Typical  Strain  Gage  (Penn  State,  2004) 


As  a  load  is  applied  the  wire  filaments,  denoted  as  etched  foil  in  Figure  20,  the 
filaments  elongate  or  shrink,  resulting  in  a  change  of  resistance  in  the  wires.  The 
resulting  difference  in  applied  voltage  to  output  voltage  can  be  used  to  calculate  the  strain 
in  the  wire. 

Resistance  Calculation  (DeLuca,  2004:  149-152) 

R  =  p  *  (L  /  A) 

Where:  p  =  Resistivety  of  the  wire 

L  =  Length  of  the  wire 
A  =  Cross  sectional  area  of  the  wire 


Calculation  of  strain  (Parga,  2004:  71) 
E  =  4  *  (Vo  /  Vs)  *  (1  /  SF) 

Where:  Vo  =  Output  voltage  of  the  bridge 

Vs  =  Voltage  applied  to  the  bridge 
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SF  =  Strain  gage  factor 


Applying  Hooke’s  Law  allows  the  stress  applied  to  be  calculated: 

£=E  *  s 

Where:  E  =  Modulus  of  elasticity  of  the  material 

Forces  and  moments  can  then  be  calculated  using: 

Force  =  a  *  s 
Moment  =  F  *  L 

When  taking  data  the  balance  will  output  six  voltages  that  are  then  interpreted  as 
forces  by  the  software.  The  forces  these  voltages  represent  are  presented  in  Table  5.  The 
maximum  load  the  balance  can  handle  in  each  direction  is  also  given.  An  understanding 
of  these  maximums  is  important  for  calibration  which  will  be  discussed  in  further 
paragraphs. 


Table  5.  Load  Limits  for  AFIT-1  Balance 


Component 

Maximum  Load 

N1  -  Normal  Force 

10  lbs. 

N2  -  Pitch  Moment 

10  in.  lbs. 

A1  -  Axial  Force 

5  lbs. 

SI  -  Side  Force 

5  lbs. 

S2  -  Yaw  Moment 

5  in.  lbs. 

LI  -  Roll  Moment 

4  in  lbs. 

Experimental  Procedure 

Before  starting  the  experiments  for  this  study  the  AFIT-1  balance  was  calibrated 
and  installed  into  the  tunnel.  Calibration  and  installation  was  done  by  Mr.  Dwight 
Gehring,  the  wind  tunnel  technician.  Steps  used  in  the  balance  calibration  are  in 
Appendix  B. 

Before  starting  the  experiments  for  this  study,  the  AFIT-1  balance  was  calibrated 
and  installed  into  the  tunnel.  Calibration  and  installation  was  done  by  Mr.  Dwight 
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Gehring,  the  wind  tunnel  technician.  Steps  used  in  the  balance  calibration  are  in 
Appendix  B. 

Once  calibrated  the  balance  and  holder  were  passed  up  through  the  bottom  of  the 
tunnel  and  attached  to  the  fixed  part  of  the  sting.  The  model  was  then  attached  to  the 
sting  with  a  mounting  block.  The  mounting  block  is  described  in  greater  detail  on  page 
74  of  the  Parga  thesis.  The  only  difference  is  that  the  mounting  block  used  for  this  test, 
Figure  22,  allowed  the  balance  to  come  into  the  block  at  a  5°  angle  relative  to  the  MAV. 
This  was  done  to  facilitate  larger  positive  elevator  deflections  of  the  tail.  Figure  21 
depicts  the  original  mounting  block  used  for  previous  studies  in  which  the  balance 
attached  to  the  block  at  an  angle  of  0°. 


Figure  21.  Mounting  Block,  Previous  Studies 


Figure  22.  Mounting  Block,  Current  Study,  5°  offset 
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Tunnel  tests  were  conducted  at  20  and  30  miles  per  hour.  The  20  mph  runs  were 
done  only  for  spoiler  characterization  to  ensure  that  the  balance  load  limits  were  not 
exceeded.  Thirty  mile-per-hour  runs  were  chosen  because  that  is  the  speed  that  was 
chosen  in  Parga’s  thesis,  thus  allowing  a  more  direct  comparison  of  results. 

To  acquire  data  a  tare  fde  was  first  created  by  running  the  model  through  a  sweep 
of  attack  angles  with  the  tunnel  motor  off.  This  gave  a  baseline  for  how  model  weight 
would  interact  with  the  various  forces  being  measured  by  the  balance.  Once  this  file  was 
saved  a  new  file  was  created  ready  to  accept  data  from  the  run  with  the  wind  on.  With 
the  fan  on  and  tunnel  up  to  speed,  the  wind  produces  the  forces  and  moments  on  the 
model  that  are  sensed  by  the  balance.  These  voltage  signals  are  then  sent  through  a  16  bit 
electronic  data  acquisition  card,  and  then  the  analog  signal  is  transfonned  to  a  digital 
signal.  Signals  are  then  amplified  and  conditioned  by  a  low  pass  filter  before  being 
stored  on  a  Pentium  computer  (Parga,  2004:  76). 

Data  for  each  flight  condition  during  the  run  was  taken  in  30  second  intervals, 
providing  enough  data  to  account  for  any  variations  due  to  vibrations  when  conditions 
were  changed.  For  the  spoiler  characterization  runs,  data  was  taken  with  the  tail  at  0° 
rotation  angle  and  with  the  tail  deflected  in  various  elevator  positions  of  -20°  to  -85°  with 
the  model  at  -  4°,  0°,  and  +  4°  angles  attack.  The  data  collection  runs  labeled  as  matrix 
runs  combined  rotation  angle  and  deflection  which  were  taken  at  +4  degrees  angle  of 
attack  and  elevator  deflection  ranging  from  +18°  to  -30°.  Rotation  angle  was  varied  from 
as  much  as  -30°  to  +30°.  Beta,  negative  yaw  angle,  runs  comprised  the  last  of  the  runs  for 
this  thesis.  Here  beta  was  swept  from  -8°  to  +8°.  The  tail  was  rotated  at  0°  and  as  far  as 
67°,  with  angle  of  attack  being  held  at  +4°. 
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Data  from  all  the  runs  was  stored  in  a  text  file  by  the  LabView  software.  These 
files  were  then  “cleaned”  by  erasing  header  and  superfluous  data  that  occurred  while  the 
tunnel  was  brought  up  to  speed.  The  clean  tab-delimited  file  was  then  called  by  a 
MATLAB  program.  The  program  used  was  initially  developed  by  Capt.  Deluca  and  Lt. 
Gebbie  (Gebbie,  2005,  1 13).  The  code  was  then  adapted  by  Parga  for  use  with  his  thesis 
(Parga,  2004:  188).  ENS  Leveron  then  readapted  it  to  include  further  revisions  to  the 
original  program  by  Lt.  Gebbie,  and  to  account  for  the  balance  being  at  a  5  degree  angle 
to  the  MAV. 

The  MATLAB®  program  takes  each  test  condition  within  the  file,  for  most  cases 
3  different  a’s,  and  averages  the  [Uco,  a,  (3,  Ni,  N2,  Si,  S2,  Ai,  t  ]  data  for  that  condition. 

This  data  is  then  used  to  calculate  the  aerodynamic  properties.  Results  are  then  written 
by  the  MATLAB®  program  to  a  file  that  is  easily  readable  by  Microsoft  Excel®.  The 
specific  equations  and  processes  used  by  the  MATLAB®  program  may  be  seen  in 
Appendix  D,  and  pages  36-45  of  DeLuca,  as  well  as  pages  77-88  of  Parga. 
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IV.  Results  and  Analysis 


Introduction 

This  chapter  presents  the  data  and  results  of  the  tests  conducted  in  the  wind  tunnel 
on  the  two  tails  discussed  in  Chapter  III.  For  data  presented  from  this  point,  tail 
deflection  as  an  elevator  is  denoted  as  5e,  and  rotation  deflection  is  denoted  as  5m.  Both 
5e  and  5m  were  measured  before  each  run  with  a  hand-held  inclinometer.  Various  trim 
positions  on  the  controller  were  noted  for  corresponding  angle  deflections  of  the  tail. 

This  allowed  multiple  5m  angles  to  be  tested  for  a  specific  5e  without  having  to  stop  the 
tunnel  and  re-measure  the  angle  for  each  data  point.  The  sign  convention  used  is  the 
same  as  that  used  in  Parga’s  thesis.  Figure  23  and  Figure  24,  illustrate  the  sign 
convention. 


Figure  23.  Elevator  Deflection  Convention 
(Parga,  2004:  59) 


Figure  24.  Rudder  Deflection  Condition 
(Parga,  2004:  59) 
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For  all  runs  the  horizontal  distance  from  the  balance  center  of  gravity  to  the  model 
center  of  gravity,  Xcm,  was  0.71  inches.  While  this  study  did  not  consider  the  movement 
of  the  center  of  gravity  of  the  model  and  its  affects  on  stability,  Lt.  Parga  did  conduct 
such  a  study.  His  results  will  be  brought  into  this  study  later.  A  pictorial  representation 
of  the  definition  of  Xcm  in  this  thesis  is  displayed  in  Figure  25. 


Figure  25.  CG  Notation  (Parga,  2004:  90) 


The  tests  run  at  20  and  30  miles  per  hour  had  the  following  average  aerodynamic 
properties: 


Table  6.  Summary  of  Flight  Conditions 


U oo  (mph) 

Mach  # 

qc(lbf  /  ft2) 

Rec 

20 

0.027 

1.09 

9.3  X  104 

30 

0.039 

2.2 

1.33  X  105 
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Additionally  roll  angle,  cp,  was  set  to  zero  for  all  test  and  therefore  P  =  -  vj/. 

The  data  is  presented  by  breaking  it  into  three  main  parts: 

a.  Characterization  of  tails  as  a  spoiler. 

b.  Characterization  of  tails  as  an  elevator  and  rudder. 

c.  Determination  of  directional  stability  of  tail  and  stabilizer  configurations. 

Spoiler  Characterization 

The  goal  of  these  experiments  was  to  detennine  how  effectively  the  tail  could  be 
used  to  rapidly  descend  in  altitude  and  act  as  an  air  brake  during  or  immediately  prior  to 
touch  down  of  the  MAV.  The  goal  of  a  spoiler  is  to  reduce  lift  and  increase  drag. 
Though  reducing  lift  will  reduce  induced  drag,  the  spoiler  should  be  capable  of 
overcoming  this  loss  with  additional  form  drag. 

During  testing  the  elevator  deflection  angle,  5e,  was  noted  to  have  some  variance 
in  it  once  the  tunnel  was  brought  up  to  speed.  This  was  due  to  tolerances  in  the 
mechanism  being  too  great  and  in  the  lack  of  torque  of  the  servo  actuators.  By  visual 
estimation,  the  elevator  deflection  angle  did  not  appear  to  become  more  positive  by  any 
more  than  5°  -  8°.  Once  the  tunnel  was  in  the  process  of  slowing  to  a  stop,  the  tail  could 
again  visually  be  seen  moving  back  to  its  pre-tunnel  run-up  position.  The  result  of  this 
action  may  be  evident  in  the  data  as  different  elevator  deflections  may  have  similar 
coefficients  or  overlap.  Figures  presented  with  respect  to  elevator  deflection  (5e)  are 
based  upon  commanded  deflection,  and  do  not  take  into  account  loss  of  deflection  angle 
due  to  the  previously  described  factors. 

Lt.  Parga  summarized  the  characterization  of  lift  with  respect  to  elevator 
deflection  angle  on  page  126  of  his  thesis.  Those  observations  are: 
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1 .  Lift  is  a  function  of  5e. 

2.  To  increase  lift  a  positive  5e  should  be  used,  similar  to  increasing  the  angle  of 
attack  of  a  wing.  To  decrease  lift  a  negative  5e  should  be  used. 

3.  Tail  rotation,  8m,  within  the  range  of  ±  20°  will  not  greatly  affect  lift 
coefficient.  Outside  of  this  range  affects  start  to  be  seen  in  the  form  of 
decreasing  lift  coefficient. 

Tail  1  was  first  run  in  the  wind  tunnel.  The  internal  linkage  allowed  this  tail  to 
have  a  maximum  deflection  of  -64°.  The  tunnel  was  first  run  at  20  mph  to  avoid 
overloading  the  balance. 


CL  due  to  Elevator  Deflection  (5e); 
20  mph,  Tail  1 


Figure  26.  CL  of  Tail  1  in  Spoiler  Configuration 

The  maximum  lift  coefficient  displayed  in  Figure  26,  which  is  for  Tail  1  at  20 
miles  per  hour  without  stabilizers,  is  1.45.  Over  all  the  tests  conducted  the  maximum  lift 
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coefficient  developed  was  1.54.  The  figure  above  starts  below  this  maximum  and 
continues  with  the  downward  trend  in  lift  that  was  summarized  in  the  Parga  thesis. 

Drag  increased  as  expected  with  a  spoiler.  The  data  did  show  the  lift  coefficient 
leveled  out  at  an  elevator  deflection  above  45°  while  the  drag  coefficient  decreased.  One 
possible  explanation  for  this  decrease  is  that  at  the  higher  deflection  angle  the  tail  was 
more  aeroelastically  defonned  more  by  the  wind  force  on  it  than  it  was  at  other  deflection 
angles.  Higher  tunnel  speeds  typically  produced  a  higher  deflection  value. 

Above  deflection  of  45°  the  Cd  takes  on  roughly  a  sine  curve  shape,  seen  in 
Figure  27.  This  is  similar  to  what  is  seen  on  a  flat  plate  at  incidence.  After 
approximately  45°  angle  of  attack  the  drag  normal  becomes  dependent  upon  sine  of  the 
deflection  angle.  A  similar  anomaly  is  found  in  the  Tail  3  at  30  mph  case.  Coefficient  of 
normal  force  on  a  flat  plate  is  further  described  to  accompany  those  figures. 


CD  due  to  Elevator  Deflection  (5e); 
20  mph,  Tail  1 


Figure  27.  CD  of  Tail  1  in  Spoiler  Configuration,  20  mph 
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The  30  mph  run  for  Tail  1  did  not  show  same  trend  in  the  decrease  of  drag 
coefficient,  Figure  28.  This  would  seem  contrary  to  the  previous  statement  of  higher 
speeds  leading  to  larger  tail  deflections.  The  tail  may  have  deflected  roughly  equally  for 
the  30  mph  runs  such  that  a  loss  in  coefficients  would  not  be  witnessed.  Lift  coefficient, 
presented  in  Figure  29,  for  the  30  mph  case  did  not  level  off  as  it  did  in  the  20  mph  case. 


CD  due  to  Elevator  Deflection  (5e); 
30  mph,  Tail  1 


Figure  28.  CD  of  Tail  1  in  Spoiler  Configuration,  30  mph 
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CL  due  to  Elevator  Deflection  (5e); 
30  mph,  Tail  1 


Figure  29.  CL  of  Tail  1  in  Spoiler  Configuration,  30  mph 

Tail  3,  with  its  different  attachment  orientation,  is  capable  of  reaching  a  maximum 
5e  of  -85°.  Thus  trends  in  the  aerodynamic  forces  and  moments  may  be  further 
characterized  beyond  what  they  were  for  Tail  1.  These  trends  continued  as  would  be 
expected  with  the  tail  projecting  more  of  its  area  perpendicular  to  the  free  stream  flow. 
However,  with  Tail  3  being  centered  vertically  a  loss  of  area  projected  into  the  free 
stream  flow  occurs.  The  MAV  fuselage  blocks  a  portion  of  this  area  where  with  Tail  1, 
being  attached  at  the  trailing  edge  of  the  wing,  does  not  experience  this  blockage  when 
deflected  in  the  negative  elevator  deflection.  This  affect  of  this  blockage  is  seen  when 
comparing  Tail  1  and  Tail  3  figures  and  noting  that  Tail  3  curves  are  shallower. 
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CL  due  to  Elevator  Deflection  (5e); 
20  mph,  Tail  3 


0  10  20  30  40  50  60  70  80  90 

Negative  Elevator  Deflection,  5e,  (degrees) 


Figure  30.  CL  of  Tail  3  in  Spoiler  Configuration,  20  mph 


CD  due  to  Elevator  Deflection  (5e); 
20  mph,  Tail  3 


Figure  31.  CD  of  Tail  3  in  Spoiler  Configuration,  20  mph 


38 


In  Figure  30  and  Figure  31,  above,  a  divergence  from  the  linear  trend  can  be  seen 
at  5e,  of  55°,  for  angles  of  attack  of  0°  and  4°.  Even  though  the  highest  lift  coefficient 
does  not  occur  at  55°  the  largest  pitch  moment  coefficient  does.  While  the  lift  coefficient 
plots  show  a  higher  lift  for  larger  alphas  the  same  is  not  true  for  the  pitch  moment. 
Blockage  by  the  model  on  the  airflow  reaching  the  tail  is  one  explanation  for  this  result. 
CG  effects  should  not  be  present  when  comparing  different  alphas  since  a  tare  file  was 
taken  to  factor  out  these  interactions. 


Cm  due  to  Elevator  Deflection  (5e); 
20  mph,  Tail  3 


Figure  32.  Cm  for  Tail  3  as  a  Spoiler,  20  mph 

The  anomaly  at  5e  of  55°  is  somewhat  similar  to  that  of  Tail  1  in  the  20  mph  case. 
First  the  deflection  of  55°  may  be  closer  to  45°  when  taking  into  account  the  deformation 
due  to  previously  described  factors.  Flat  plate  tests  conducted  by  Sighard  Hoerner  found 
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the  coefficient  of  normal  force  was  found  to  increase  linearly  until  approximately  45°. 
After  this  angle  a  sharp  drop  of  the  normal  force  coefficient  to  a  constant  value  is  noticed. 
Above  45°,  Cn0rmai  becomes  constant  with  Cl  normal  becoming  a  function  of  cosine  a  and 
Cd  normal  a  function  of  sine  a. 

Cntrml 


Figure  33.  Normal  Force  Coefficient  on  a  Square  Plate  (Floerner,  1965:  3-16) 

This  reference  does  prove  beneficial  in  determining  the  anomalies  found  in  Figure 
30,  Figure  31,  and  Figure  32.  The  equations  for  Cl  normal  and  Cd  normal  show  that 
beyond  45°  Cl  would  continue  to  decrease  to  0  at  90°  while  Cd  would  increase,  reaching 
its  maximum  value  at  90°.  The  fact  that  Cn0nmi  reaches  a  maximum  then  drops  fits  the 
shape  of  the  Cm  plot,  Figure  32.  The  largest  moment  would  be  expected  to  correspond 
with  the  largest  forces  on  the  tail.  Additionally  the  pitch  moment  arm  changes  with 
respect  to  the  elevator  deflection  angle.  For  a  constant  normal  force  applied  to  the  tail  the 
resulting  moment  about  the  CG  would  decrease  as  deflection  increases  through  45°.  This 
would  give  a  similar  shape  to  the  pitch  moment  curves  as  depicted  in  Figure  32. 
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Another  interesting  result  is  that  the  similar  anomalies  were  not  found  at  30  mph. 
In  this  case,  peak  value  of  Cm  was  found  at  65°.  Cl  showed  the  greatest  change  from  the 
previous  deflection  value  while  CD  showed  a  similar  trend.  However  where  Cl  leveled 
off,  Cd  continued  to  rise;  again  this  holds  fairly  consistent  with  the  equations  of  each  as  a 
trigonometric  function  of  alpha. 


CL  due  to  Elevator  Deflection  (5e); 
30  mph,  Tail  3 


Figure  34.  CL  for  Tail  3  as  a  Spoiler,  30  mph 
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CD  due  to  Elevator  Deflection  (5e); 
30  mph,  Tail  3 


Figure  35.  CD  of  Tail  3  as  a  Spoiler,  30  mph 


Cm  due  to  Elevator  Deflection  (5e); 
30  mph,  Tail  3 


Figure  36.  Cm  of  Tail  3  as  a  Spoiler,  30  mph 
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Overall  the  two  tails  gave  similar  results  for  the  spoiler  tests.  The  added  increase 
in  deflection  was  needed  to  lower  the  lift  coefficient  to  the  value  that  Tail  1  achieved  at  a 
less  negative  angle.  The  general  shapes  of  the  coefficient  plots  varied  while  their  values 
stayed  roughly  similar.  Cl  and  Cd  of  the  two  tails  generally  did  not  differ  by  any  more 
than  0.1,  while  Cm  varied  by  0.2  or  slightly  more. 

The  attachment  of  the  two  tails  no  doubt  had  an  effect  on  the  pitching  moment. 
Referring  to  the  attachment  differences  in  the  tails  described  in  Chapter  III,  Tail  3 
experienced  more  blockage  of  its  area  from  the  fuselage  than  Tail  1.  As  discussed,  the 
attachment  also  reduced  the  length  of  the  tail,  thus  possibly  accounting  for  some  of  the 
loss  of  moment  that  is  evident  in  Figure  37.  Both  tails  did  function  as  spoilers  and  would 
be  suitable  to  act  as  an  air  brake  upon  landing.  Tail  3  gave  the  highest  drag  for  a  large 
negative  deflection  angle. 
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Tail  1  vsTail  3;  CL,  CD  and  Cm 
20  mph,  alpha  =  4 


Figure  37.  Comparison  of  Tail  1  and  3  Selected  Properties 

Rudder  and  Elevator  Characterization 

All  wind  tunnel  tests  presented  in  the  current  section  were  conducted  at  an  angle 
of  attack  of  4°,  P  =  0°,  and  at  30  mph.  Similar  test  were  conducted  by  Lt.  Parga. 
However,  tests  here  sought  to  extend  the  5m  range.  His  test  were  conducted  at  and 
between  5m=  +/-  20°,  while  rotation  angles  in  this  study  were  up  to  +/-  30°. 

Data  was  collected  in  a  matrix  form  similar  to  Table  7  and  Table  8.  Again 
rotation  and  deflection  angles  were  measured  with  a  hand-held  inclinometer  prior  to  the 
test  and  angles  were  matched  with  specific  trim  positions  on  the  hand-held  controller. 
Differences  in  the  rotation  angles  used  were  a  result  of  the  tail  being  used.  Each  tail  was 
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capable  of  different  rotation  angles  based  upon  its  attachment  and  constant  tweaking  of 
the  internal  control  linkages  while  changing  configurations  between  tests. 

Table  7.  Test  Matrix  for  Tail  3 


5m 


Table  8.  Test  Matrix  for  Tail  1 

5m 


Additionally,  these  tests  were  performed  to  determine  whether  the  rotatable  tail  is 
a  worthy  substitute  for  a  rudder  and  is  capable  of  supplying  similar  forces  and  moments 
as  a  conventional  rudder.  The  range  of  rotation  angles,  for  Tail  3,  fall  within  the  range  of 
most  conventional  rudders  as  they  are  often  limited  to  +/-  30°  because  beyond  this  range 
effectiveness  falls  off  abruptly  (Perkins,  1958:  329). 

Tail  effects  on  Lift  coefficient 

Results  are  similar  to  those  in  the  spoiler  characterization  section  and  to  those  in 
the  thesis  of  Parga,  page  126.  Figure  38  and  Figure  39  show  lift  coefficient  to  be  a 
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function  of  elevator  deflection.  It  is  not  affected  greatly  by  rotation  angles  inside  of  -  20° 
<  5m  <  +  20°.  Above  the  absolute  value  of  20°,  a  slight  drop  off  in  lift  starts  to  be 
noticeable.  As  the  tail  rotates  the  projected  horizontal  area  is  becoming  a  vertical 
projected  area.  With  area  that  once  was  producing  lift  now  producing  a  side  force,  loss  of 
lift  is  expected. 

Negative  elevator  deflection  gives  a  decrease  in  lift  coefficient  from  the  neutral 
position  Cl,  while  positive  deflection  gives  an  increase  in  lift  coefficient.  This  is 
consistent  with  a  wing  or  flat  plat  producing  higher  lift  when  at  a  positive  angle  of  attack 
and  less  lift  when  at  a  negative  angle  of  attack. 

All  contour  plots  are  for  a=  4°;  which  was  chosen,  not  because  of  maximum  L/D, 


but,  because  the  work  of  Deluca  and  Parga  was  performed  at  this  angle  of  attack.  This 


allowed  a  comparison  to  be  easily  drawn  between  values  obtained  from  one  study  to  the 
next. 
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Figure  38.  CL  Contour  Plot,  Tail  1 
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CL  Contour  Plot 
Tail  3 


Figure  39.  CL  Contour  Plot,  Tail  3 

Tail  Effects  on  Drag  Coefficient 

As  would  be  expected  and  presented  in  the  spoiler  section,  drag  coefficient  is  a 
function  of  elevator  deflection  angle.  The  Tail  1  plot,  Figure  40,  shows  no  dependence  of 
drag  on  8m.  Plots  of  both  tails  show  a  “drag  bucket”  where,  as  elevator  deflection  is 
increased  from  either  side  of  zero,  the  increase  in  parasitic  drag  overcomes  the  loss  of 
induced  drag  for  a  net  incase  causes  an  total  drag. 

A  negative  5e  results  in  loss  of  lift,  as  witnessed  in  Figure  38  and  Figure  39.  This 
loss  of  lift  also  results  in  loss  of  induced  drag,  however;  drag  is  seen  to  increase  as  5e 
increases.  The  tail  acting  as  a  spoiler  and  increasing  form  drag  is  attributed.  For  a 
positive  5e  lift  is  seen  to  increase  once  5e  drops  below  ~10°.  The  corresponding  increase 
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in  induced  drag  is  attributed  to  the  increase  seen  in  Figure  40  and  Figure  41,  for 


increasing  positive  5e. 
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Figure  40.  CD  Contour  Plot,  Tail  1 
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CD  Contour  Plot 
Tail  3 
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Figure  41.  CD  Contour  Plot,  Tail  3 

Tail  Effects  on  Pitch  Moment  Coefficient 

Similar  to  the  characterization  in  Parga’s  thesis  on  page  132,  pitch 
moment  was  again  found  to  be  dependent  upon  elevator  deflection  with  little  influence 
due  to  5m.  Positive  elevator  deflection  gave  a  negative  increase  in  moment  coefficient. 
Similarly,  a  negative  elevator  deflection  gives  a  positive  increase  in  pitch  moment 
coefficient.  Values  for  the  moment  coefficient  are  different  from  those  of  Lt.  Parga. 
Figure  99  of  Parga’s  thesis,  page  133,  shows  Cm  values  in  the  range  of  0.0175  to  -0.2636 
while  the  current  study  values  in  the  range  of  -0.05  to  -0.5  were  attained  with  larger 
deflection  angles.  This  is  attributed  to  a  different  location  of  the  CG  of  the  MAV  but  also 
position  of  the  tail  with  respect  to  vertical  position  of  the  downwash  from  the  wing  must 
be  brought  into  consideration.  Additionally,  the  new  balance  block  introduced  some 
affects  on  location  of  balance  CG  to  MAV  CG  that  were  ignored.  The  wing  was  also 
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removed  and  reattached.  Even  though  it  was  placed  in  its  original  position  any  difference 
in  incidence  angle  or  location  would  have  an  effect  of  differing  values  calculated  from 
this  and  the  previous  study  by  changing  downwash  from  the  wing.  However,  values 
obtained  may  be  changed  by  varying  the  CG  location.  The  range  of  values  obtained  is 
the  main  area  of  interest. 

The  data  presented  does  state  that  the  MAV  will  have  a  problem  trimming 
because  the  elevator  will  have  to  be  deflected  to  a  considerable  negative  angle  which 
would  be  detrimental  to  efficient  flight.  Based  upon  the  dependence  of  Cm  to  CG 
location,  changing  the  CG  would  be  the  best  course  of  action  to  achieve  desirable 
moment  coefficients. 
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Figure  42.  Cm  Contour  Plot,  Tail  1 
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Figure  43.  Cm  Contour  Plot,  Tail  3 


Tail  Effects  on  Yaw  Moment  Coefficient 

As  mentioned  earlier  in  this  chapter,  the  previous  study  on  this  subject  did  not 
show  the  rotatable  tail  measuring  up  to  traditional  rudders  in  yaw  moment  coefficient. 

Lt.  Parga  showed  a  Cn  range  of  -0.013  to  +0.013  for  tail  1.  A  typical  rudder  has  a  Cn 
range  of  -0.02  to  +  0.02  (Barlow  et.  ah,  1999:  527). 

Behavior  exhibited  by  the  tail  on  Cn,  is  again  similar  to  previous  observations. 
(Parga;  2004).  Tests  conduced  pointed  to  the  following  behavior  of  Cn  with  respect  to 
tail  movements. 

1 .  Cn  is  a  combination  of  5e  and  Sm. 

2.  Positive  5e  and  5m  is  seen  to  produce  more  favorable  Cn  values  for  a  turn.  This 
is  similar  to  some  observations  of  birds  using  a  depressed  and  twisted  tail  to 
execute  a  level  turn. 
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3.  Control  reversal  is  possible.  For  a  given  5e  and  8m  a  value  of  Cn  is  achieved. 


Changing  5e  only  will  result  in  a  change  in  the  sign  of  the  Cn  value.  Thus 


leading  to  a  different  yaw  than  was  originally  commanded. 
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Figure  44.  C„  Contour  Plot,  Tail  1 
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Figure  45.  C„  Contour  Plot,  Tail  3 

The  greatest  change  in  Cn  range  is  seen  in  Figure  45  for  Tail  3,  as  it  had  the 
highest  rotation  angles,  -30°  to  +30°.  Tail  3  was  capable  or  reaching  a  Cn  value  of  -0.02 


however  it  reached  a  slightly  less  positive  value  0.016.  This  showed  that  the  rotatable  tail 


is  capable  of  delivering  essentially  the  same  values  as  a  typical  rudder.  Extending  the 


range  of  Cn  is  possible,  and  will  be  discussed  in  a  future  section  on  improving  stability. 


Figure  46  represents  a  rudder  power  curve  for  a  traditional  rudder.  The  data  from 


this  figure  is  plotted  in  Figure  47  and  Figure  48  against  values  obtained  for  Tail  1  and 


Tail  3.  The  severity  of  the  elevator  deflection  greatly  affects  the  rudder  power  the  tail 


generates.  An  equivalent  range  and  slope  to  that  seen  Figure  46  is  achieved  by  the 


rotatable  tail  at  the  large  negative  elevator  deflection. 
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Figure  46.  Rudder  Power  Curve  (Perkins,  1958:  330) 


Rudder  Power,  Tail  1 


Figure  47.  Rudder  Power  Curve  for  Tail  1 
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Rudder  Power,  Tail  3 


Figure  48.  Rudder  Power  Curve  for  Tail  3 


Tail  Effects  on  Roll  Moment  Coefficient 

Similar  characteristics  were  found  between  prior  results  and  data  taken  for  this 
test  (Parga;  2004).  The  rotatable  tail  did  not  display  the  amount  of  roll  control  that  is 
present  in  aileron  control  surfaces.  Tail  3  does  not  display  any  noteworthy  differences 
over  Tail  1  as  far  as  roll  moment  coefficient  is  concerned.  The  effects  of  both  rotatable 
tails  are  characterized  by  roll  moment  being  a  function  of  both  5e  and  5mwith  the  largest 
achieved  values  being  at  the  extremes  of  each  commanded  angle,  as  evidenced  in  Figure 
49  and  Figure  50. 
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Figure  49.  Ci  Contour  Plot,  Tail  1 
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Figure  50.  Q  Contour  Plot,  Tail  3 


Directional  Stability  Investigation 

As  presented  earlier,  the  aircraft  configured  with  the  rotatable  tail  exhibits 
problems  with  directional  stability.  An  attempt  to  alleviate  this  problem  was  carried  out 
in  research  done  by  Parga  by  attaching  vertical  stabilizers  to  a  rotatable  tail.  This  was 
found  to  improve  stability  some  but  not  enough  to  be  considered  stable.  Reviewing  the 
stabilizer  used  in  this  study  it  consisted  of  an  area  of  7. 13  in2  attached  at  the  lower  aft  end 
of  the  fuselage.  The  vertical  tail  coefficient  for  the  stabilizer  is  0.04  for  the  tail  with  8e 
and  Sra  =  0°,  and  the  horizontal  tail  coefficient  is  0.34.  The  tail  coefficient  values  would 
change  with  a  change  in  8m.  If  8m  was  deflected  to  90°,  horizontal  tail  coefficient  would 
be  0  while  vertical  tail  coefficient  would  be  greater  than  .34,  due  to  the  contribution  of 
tail  area  being  projected  vertically.  Chapter  3  describes  the  stabilizer  more  completely 
and  Table  2  gives  stabilizer  properties. 

Data  was  taken  by  varying  P  angles.  P  had  the  range  of  -8°  to  +8°  for  this  test. 
Angle  of  attack  was  held  constant  at  4°  and  the  tunnel  speed  was  30  mph.  Data  was  taken 
for  both  Tail  1  and  Tail  3,  each  with  and  without  the  stabilizer  attached.  Some  of  the 
following  plots  are  duplicates  of  those  presented  earlier;  this  is  done  for  ease  of 
comparison.  The  stabilizer  had  little  effect  on  the  tail  Cl,  Cm,  and  Cd-  The  stabilizer  no 
doubt  added  drag,  but  did  not  affect  the  tail  contribution  of  drag  other  than  causing  minor 
flow  interference  at  the  positive  angles  of  elevator  deflection. 
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Figure  51.  Cn  Contour  Plot,  Tail  1 


Figure  52.  Cn  Contour  Plot,  Tail  1 
w/stabilizer  s 


Tail  1  Comparision  with  and  with  out  Stabilizer 
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Figure  53.  Tail  1  Stabilizer  Comparison 
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The  stabilizer  did  have  the  effect  of  increasing  the  range  of  Cn.  The  properties  of 
the  tail  with  respect  to  Cn  were  kept  the  same  as  without  a  stabilizer.  The  stabilizer’s 
affect  on  Tail  3  again  served  to  increase  the  Cn  range.  It  did  however  slightly  decrease 
the  positive  range  from  0.016  to  0.01.  Cn  properties  remained  the  same  as  discussed  in 
the  elevator/rudder  section.  The  stabilizers  did  have  the  affect  of  allowing  a  larger  value 
of  yaw  coefficient  to  be  reached  for  each  incremental  increase  in  6e  and  5m.  Both  plots, 


Figure  54  and  Figure  55,  show  how  elevator  and  rotation  deflection  must  be  used 


together  to  achieve  an  effective  yaw  moment. 
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Figure  54.  Cn  Contour  Plot,  Tail  3 


Figure  55.  Cn  Contour  Plot,  Tail  3 
w/stabilizers 


Comparison  charts  show  that  it  is  possible  for  both  tails  to  give  the  same  yaw 
moment  coefficient  values  as  a  traditional  rudder  despite  the  relative  small  tail  volume 
coefficients.  A  typical  rudder  will  provide  a  Cn  of  ±  0.02.  The  configurations  and 
deflections  used  in  this  study  produced  maximum  absolute  values  of  0.018  and  0.03. 
Figure  56  allows  a  comparison  between  a  typical  rudder  and  the  tails  tested  to  be  drawn. 
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The  range  and  attainable  values  are  shown  to  be  similar.  While  the  slope  of  Cn  for  the 
tested  tails  is  not  as  great  as  that  of  a  typical  rudder,  the  attainable  values  are  of  more 
interest.  Even  though  an  additional  10°  of  rotation  is  needed  the  desirable  value  of  0.02  is 
reached. 


Tail  3  Comparison  with  and  with  out  Stabilizer 


Figure  56.  Tail  3  Stabilizer  comparison 


The  stabilizer  had  a  mixed  effect  with  regards  to  directional  stability  at  a  cruise 
condition,  5e=  0°.  The  stabilizer  did  have  the  affect  of  improving  stability  for  both  Tail  1 
and  Tail  3  when  no  rotation  deflection  is  present.  When  a  5m  is  present  the  stabilizers  had 
the  affect  reducing  the  Cn  maximum  values.  However,  the  slope  of  the  Cn  curve  was 
increased.  The  stabilizers  generally  shifted  the  Cn  curve  down  to  what  a  typical  curve 
would  look  like  with  0°  yaw  producing  a  zero  Cn  value.  The  data  presented  in  Figure  57 
and  Figure  58  show  that  with  and  without  stabilizer,  the  MAY  would  be  more 
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directionally  stable  in  a  turn,  i.e.  with  a  5m  of  a  value  other  than  zero.  A  typical  data  set 
for  both  directional  and  roll  stability  is  shown  in  Figure  59. 


Directional  Stability  for  Tail  1  Configuratons 
5e  =  0 


Figure  57.  Directional  Stability  Tail  1  Configurations 
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Directional  Stability  for  Tail  3  Configurations 
5e  =  0 


Figure  58.  Directional  Stability  Tail  3  Configurations 
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Figure  59.  Typical  Yaw  Characteristics  (Barlow  et  al,  1999:  531) 
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Roll  Stability  of  Tail  1  Configurations 


Roll  Stability  for  Tail  3  Configurations 
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Figure  64.  Q  Contour  Plot,  Tail  3 


Figure  62.  Ci  Contour  Plot,  Tail  1 
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Cx  Contour  Plot 
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Figure  63.  Ci  Contour  Plot,  Tail  1 
w/stabilizers 


Figure  65.  C,  Contour  Plot,  Tail  3 
w/stabilizers 


The  roll  moment  coefficient  is  only  slightly  affected  by  the  rotation  of  the  tail. 
Common  values  are  within  the  range  of  -0.02  to  +  0.02,  when  an  aileron  is  used  (Barlow 
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et  al,  1999:  500).  The  range  of  values  obtained  with  the  rotatable  tail  was  about  ten  times 
smaller,  as  shown  in  Figures  60-65. 

Stabilizers  had  little  effect  on  roll  moment  coefficient  properties  and  more  on  the 
range  of  Q  value.  Clearly,  from  comparing  the  four  figures  directly  above,  the  effect  on 
Tail  1  is  extending  the  attainable  values.  The  small  offset  from  Q  =  0  has  also  been 
noted  by  DeLuca  and  Parga.  One  possible  reason  may  be  that  a  slight  offset  in  the  wing 
attachment  exists. 

The  directional  stability  of  the  MAV  is  probably  the  area  of  greatest  improvement 
and  offers  the  most  improvement  for  future  development.  As  stated  before  the  original 
MAV  with  the  rotatable  tail  had  directional  stability  problems.  The  work  of  Parga  added 
stabilizers  to  the  tails  to  address  this  problem,  but  this  complicated  the  results.  Adding 
stationary  vertical  stabilizers  in  the  manner  of  this  study  again  improved  stability  but  not 
enough  to  be  considered  a  satisfactory  stable  model.  The  lack  of  stability  would  require 
automatic  controls  if  it  were  to  be  flown  by  a  human  pilot.  Improving  stability  to  an 
acceptable  level  would  allow  a  human  to  pilot  the  model  with  simply  a  remote  control 
transmitter. 

Static  stability  derivatives  calculated  by  Lt.  Parga  are  displayed  in  Table  9. 


Table  9.  Static  Stability  Derivatives  by  Lt.  Parga 


Derivative 

Tail  1 

Tail  2 

Original  UAV 

Generic 

Value 

dCmlda* 

-0.0166 

-0.01 

-0.0466 

-0.0119 

d  Cn  /  <9(3 

-0.0003 

0.0007 

0.00078 

0.00123 

d  C,  /  <9(3 

-0.0018 

-0.0017 

-0.00076 

-0.00129 

*  CG  dependent 


65 


Figure  66  represents  the  correlation  of  Cn  and  Ci  stability  derivatives  with  respect  to  yaw 
and  their  affect  on  stable  flight. 
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Figure  66.  Stability  Combinations  of  C„,|,  and  Ci¥  (Barlow  et  al,  1999:  530) 


Changing  the  derivatives  in  Table  9  so  they  are  with  respect  to  yaw  angle,  will 
allow  a  more  direct  comparison  with  the  data  presented  in  this  study,  in  which  derivatives 
are  given  with  respect  to  yaw  angle.  Similar  stability  derivates  found  in  this  test  for  both 
the  stabilizer  on  and  off  conditions  are  presented  in  the  following  table,  Table  10,  for 
which  elevator  and  rotation  deflection  were  zero. 


Table  10.  Static  Stability  Derivatives 


Derivative 

Tail  1 

Tail  3 

Tail  1  w/stabilizer 

Tail  3  w/stabilizer 

d  Cn  /  dip 

0.0002 

0.0002 

-0.0005 

-0.0005 

d  C|  /  dip 

0.0022 

0.0019 

0.0022 

0.0022 

Recalculating  the  derivatives  with  a  rotation  angle  present  allows  the  fact  that  the 
MAV  is  more  stable  in  a  turn  than  when  in  straight  level  flight  to  be  witnessed.  Table  11, 
below,  is  one  condition  for  Tail  3  with  stabilizers  and  rotation  deflection  of  54°. 
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Table  11.  Static  Stability  Derivatives  for  5,.n  =  54“ 


Derivative 

Tail  3 

3  Cn  /  di|J 

-0.0009 

3  C|  /  dijj 

0.0018 

Plotting  values  from  each  of  the  above  tables  onto  Figure  71  gives  an  idea  of  how  the 
MAY  performs  with  respect  to  directional  stability. 
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Figure  67.  Plotted  Stability  Derivatives  (Barlow  et  al,  1999:  530) 


Taking  the  above  figure  and  the  relative  “weakness”  of  the  stability 
derivatives  leads  to  one  conclusion  about  the  stability  of  the  MAV.  The  CG  of  the  MAV 
is  too  far  aft  or  rather  too  much  of  the  model  is  forward  of  the  CG.  Parga  mentioned  this 
in  his  thesis  and  conducted  a  study  on  how  CG  affects  the  MAV  stability. 
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Limitations  of  Experimental  Effort 


No  matter  how  meticulous  one  is  in  set  up  and  data  acquisition,  errors  will  always 
be  present.  Numerous  sources  of  error  in  this  study  are  easily  identifiable.  The  first  error 
source  was  discussed  in  previous  chapters  with  the  drift  in  tail  8e  due  to  dynamic  pressure 
and  its  actuation  and  structural  limits.  Additionally,  deflection  and  rotation  angle  drifts 
were  present  when  using  the  remote  controller  to  adjust  the  tail  position.  Ideally  moving 
controller  settings  from  one  position  to  another  and  then  back  again  should  have 
produced  the  original  commanded  deflection  or  rotation.  However  this  was  not  the  case 
due  to  errors  within  the  hand  held  controller,  servo  actuation  system,  and  tail  attachment 
all  compounding  upon  one  another.  Lt.  Parga  suggests  variations  of  -2°  to  +2°  in 
controller  settings.  This  test  found  that  similar  variations  were  present. 

Balance  accuracy  is  another  source  of  error.  The  AFIT-1  balance  limited 
resolution.  Reviewing  Appendix  B  and  the  balance  calibration  plots,  shows  that  the 
accuracy  of  the  balance  is  consistent  throughout  the  range  of  loadings.  Table  12 
describes  the  resolution  of  the  balance  for  the  varying  forces  and  moments  it  measures. 


Table  12.  Balance  Error 


Measurement 

Resolution  (Ibf) 

v,. 

Normal  Force 

0.012 

0.12% 

Axial  Force 

0.002 

0.04% 

Side  Force 

0.0035 

0.07% 

Pitch  Moment 

0.005 

0.05% 

Roll  Moment 

0.0044 

0.11% 

Yaw  Moment 

0.0035 

0.07% 

The  tail  attachment,  discussed  in  Chapter  III,  is  not  perfectly  symmetric.  This 
accounts  as  a  partial  explanation  for  some  of  the  forces  and  moments  not  being  zero  or 
centered  about  5e  and  5m  of  zero  degrees.  Figure  68  taken  from  Lt.  Parga  displays  the 
asymmetric  attachment  of  the  tails. 
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Figure  68.  Asymmetric  Attachment  of  Tails 


Minute  mathematical  errors  were  introduced  by  MATLAB  which  used  a  4th 
degree  polynomial  to  calculate  tare  values  (Parga,  2004:  159).  The  MATLAB  software 
had  a  floating  point  precision  2.22  X  10'16.  Using  a  4th  order  polynomial  vice  a  5th  order 
introduced  an  errors  starting  at  the  1  X  10'5  decimal  place.  Errors  introduced  into  the  test 
signal  by  the  quantization  of  the  analog  to  digital  converter  16  bit  acquisition  card  had  an 
accuracy  of  ±1.53  micro  volts  (DeLuca,  2004:  131).  These  error  sources,  though  present, 
are  insignificant. 
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V.  Conclusion 


This  study  sought  to  determine  whether  the  rotatable  tail  can:  be  used  to  aid  a 
survivable  recovery;  provide  a  yaw  moment  equal  to  that  of  a  traditional  tail 
configuration;  and  be  utilized  on  a  directionally  stable  aircraft.  Extending  the  range  of 
deflections  added  to  data  and  plots  that  were  produced  in  that  theses  and  provided  a  more 
complete  picture  of  the  characteristics  of  the  rotatable  tail.  It  was  noted  that  the  physical 
values  obtained  were  slightly  different  from  those  reported  earlier  (Parga;  2004).  A 
difference  in  the  CG  reference  position  has  been  noted  as  a  probable  cause.  Additionally, 
changes  in  wing  position  and  incidence  angle  to  the  fuselage  may  attribute  to  these 
differences.  The  range  of  obtained  values  remains  the  most  important  factor  in  tests 
conducted.  Extending  the  deflection  angles  resulted  in  an  increase  of  the  range  of 
achievable  yaw  moment  coefficient  values. 

The  research  completed  by  Parga  concluded  that  the  rotatable  tail  was  capable  of 
providing  roughly  65%  of  the  yaw  control  force  of  a  typical  rudder.  This  study  proved 
that  the  rotatable  tail  configuration  was  capable  of  providing  yaw  control  equivalent  to  a 
nonnal  rudder.  Roll  control  was  less  than  the  typical  desired  value,  but  was  comparable 
to  that  achieved  with  the  original  MAV  with  the  V-Tail  (DeLuca;  2004). 

The  rotatable  tail  with  large  negative  5e,  less  than  negative  40°,  changes  from  a 
streamlined  lifting  surface  to  a  bluff  body  producing  drag  to  slow  the  velocity  of  the 
MAV.  Research  conduced  showed  that  the  general  trend  from  increasing  negative  5e  was 
a  reduction  in  Cl  and  an  increase  in  Cd. 
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By  extending  the  range  of  elevator  and  rotation  deflection  (5e  and  5m),  the 
rotatable  tail  is  capable  of  delivering  a  yaw  similar  to  that  of  a  traditional  redder,  despite 
a  small  tail  volume  coefficient.  However,  the  coupled  controls  mandate  a  more  complex 
flight  control  system. 

The  large  area  of  the  model  forward  of  the  CG  counteracts  the  stabilizing  affects 
the  tail  and  stabilizer  produces.  A  longer  moment  arm  is  needed  between  the  tail  and  the 
CG  to  improve  directional  stability,  and  is  recommended.  The  large  elevator  deflection 
needed  to  trim  the  MAV  for  level  flight  also  confirms  this  conclusion.  To  achieve  a  more 
forward  CG,  the  wing  should  also  shift  forward. 

This  solution  would  allow  for  the  increase  in  moment  arm  desired  while  negating 
pitch  and  trim  problems  of  having  the  CG  far  forward  of  the  wing. 

Suggestions  for  future  work  include: 

1 .  Developing  a  mixing  algorithm  to  account  for  or  use  the  control  reversal 
to  an  advantage  in  a  remote  control  device. 

2.  Moving  the  wing  forward  to  allow  for  more  tail  effectiveness. 

3.  Consider,  increasing  the  size  of  the  tail,  thereby  slightly  increasing  length 
of  the  moment  arm  and  area  that  is  capable  of  providing  control. 

4.  Consider,  increasing  stabilizer  size  to  improve  stability. 
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Appendix  A:  Data  Tables 

*Tail  3  is  referred  to  as  Tail  2  in  data  tables 


Table  13.  Spoiler  Run  Data,  Tail  1  w/stabilizer 


Beta 

Delta  e 

Delta  rn 

Tail 

stabalizers 

M# 

Re# 

q_c 

Uoo 

alpha_c 

CL 

C_D_c 

Cl_cg_w 

Cm_cg_c_ 

Cn_cg_w 

C_Y\r 

0 

64 

0 

1 

1 

0.02635 

89826.4 

1.0004 

20.4124 

-4.9335 

0.15045 

0.34857 

0.00247 

0.48379 

-0.006 

0.00263 

0 

64 

0 

1 

1 

0.02632 

89730 

0 . 99825 

20.3905 

0.48235 

0.79735 

0.36569 

0.00652 

0.33931 

-0.0054 

-0.015 

0 

64 

0 

1 

1 

0.02622 

89383.7 

0 . 99056 

20.3118 

4.7732 

1.25864 

0.42961 

0.00884 

0.17556 

-0.0052 

-0.0282 

0 

20 

0 

1 

1 

0.02637 

89905.8 

1.00217 

20.4305 

-4.8147 

0.41941 

0.12595 

0.0033 

0.17964 

0.00166 

-0.0209 

0 

20 

0 

1 

1 

0.02633 

89760.6 

0 . 99893 

20.3975 

0.61074 

1.08789 

0.14669 

0.00767 

-0.0143 

0.00207 

-0.037 

0 

20 

0 

1 

1 

0.02622 

89380.4 

0 . 99049 

20.3111 

4.89874 

1.54273 

0.23084 

0.01092 

-0.2412 

0.0018 

-0.0492 

Table  17.  Matrix  Run  Data,  Tail  1 


Beta 

5e 

5m 

Tail 

stabalizers 

M# 

Re# 

q_c 

Uoo 

alpha_c 

CL 

Cd 

C_, 

Cm 

Cn 

C_Y\r 

0 

0 

23 

1 

0 

0.03925 

133798 

2.21963 

30.4058 

4.77815 

1.26986 

0.14408 

0.0059 

-0.3418 

-0 . 0019 

-0.0236 

0 

0 

14 

1 

0 

0.03927 

133896 

2.22287 

30.4279 

4.8677 

1.27562 

0.14519 

0.00598 

-0.3553 

0.00025 

-0.0291 

95!^V 

0.03924 

133785 

2.21918 

30.4027 

4.86959 

1.27989 

0.14571 

0 . 00623 

-0.3588 

0.00174 

-0.0334 

133699 

2.21635 

-0.3571 

133808 

2.21996 

1.2765 

0.1453 

-0.3538 

-0.0376 

0 

0 

-12 

1 

0 

133699 

2.21633 

-0.3474 

-0.04 

0 

-10 

23 

1 

0 

133560 

2.21173 

1.2068 

-0.2211 

-0.0516 

0 

-10 

14 

1 

0 

133475 

2.20893 

-0.0394 

0 

-10 

8 

1 

0 

133471 

2.20878 

-0.036 

0 

-10 

1 

0 

133523 

2.21052 

1.2097 

-0.2217 

133646 

2.21458 

-0.0001 

-0.0277 

0 

-0.21 

-0 . 0014 

0 

-20 

23 

1 

0 

133896 

2.22286 

-0.0781 

0 

-20 

saaBi 

133816 

2.22023 

-0.0853 

0.01115 

-0.0578 

0 

-20 

£g^^B 

133768 

2.21862 

0.00513 

-0.0467 

0 

-20 

aa^Bi] 

134137 

2.23089 

■WiIilcMi 

0 

-20 

133587 

2.21261 

0 

-20 

-12 

i 

0 

133605 

2.21322 

0 

-30 

23 

i 

0 

133688 

2.21598 

-0.0794 

0 

-30 

14 

i 

0 

133889 

2.22263 

-0.0638 

0 

-30 

8 

i 

0 

133774 

2.21884 

-0.0528 

0 

-30 

i 

0 

133777 

2.21892 

1.1115 

0 

-30 

-7 

i 

0 

133838 

2.22094 

-0.0097 

0 

-30 

133845 

2.22117 

0 

18 

23 

i 

0 

133828 

2.22061 

-0.5019 

-0 . 0161 

0 

SIGHED 

ll^B] 

fl^Bl 

133784 

2.21917 

1.36761 

-0.5209 

-0 . 006 

-0.0134 

IB^BD 

B^B] 

133820 

2.22035 

30.4107 

0.1715 

-0 . 0016 

y^BE 

IB^BS 

IH^^D 

fB^^BD 

133856 

2.22156 

30.419 

■nmutl 

-0.0361 

gBH] 

0.03925 

133798 

2.21962 

30.4057 

4.90738 

1.36541 

0.17166 

0.00848 

-0.5187 

-0.0507 

0 

18 

-12 

i 

0 

0.03926 

133859 

2.22166 

30.4196 

4.90432 

1.35848 

0.17111 

0.00911 

-0.5106 

0 . 0121 

-0.0641 

Table  18.  Matrix  Run  Data,  Tail  1  w/stabilizer 


Yaw 

5e 

5m 

Tail 

stabalizers 

M# 

Re# 

q_c 

Uoo 

alpha_c 

CL 

Cd 

c_, 

Cm 

Cn 

C_Y\r 

0 

0 

-20 

1 

1 

0.03839 

130876 

2 . 12373 

29.7416 

4.81178 

1.34594 

0.15763 

0.00853 

-0.3593 

0.00299 

-0.0471 

0 

0 

-12 

1 

1 

0.03839 

130897 

2 . 12442 

29.7465 

4.89991 

1.3485 

0.15803 

0.00851 

-0.3605 

0.00261 

-0.0462 

0 

0 

0 

1 

1 

0.03841 

130941 

2 . 12584 

29.7564 

4.90119 

1.3514 

0.15734 

0.00837 

-0.3649 

0.00203 

-0.0453 

0 

0 

7 

1 

1 

0.0384 

130920 

2 . 12517 

29.7517 

4.89991 

1.34851 

0.15684 

0.00826 

-0.3596 

0 . 0014 

-0.0438 

0 

0 

14 

1 

1 

0.0384 

130928 

2 . 12543 

29.7536 

4.89868 

1.34573 

0.15649 

0.00821 

-0.3546 

0.00115 

-0.0435 

0 

0 

25 

1 

1 

0.03838 

130855 

2 . 12306 

29.737 

4.89758 

1.34323 

0.15642 

0.0083 

-0.3489 

0 . 001 

-0.0434 

0 

-10 

-20 

1 

1 

0.03834 

130721 

2 . 11872 

29.7066 

4.78136 

1.27712 

0.16489 

0.00621 

-0.2273 

-0 . 0048 

-0.025 

0 

-10 

-12 

1 

1 

0.03837 

130822 

29.7293 

4.8678 

1.27584 

0.16517 

0.00647 

V.  MM 

-0 . 0029 

-0.0302 

iSflfiE?] 

IS^^K] 

ic^b 

BC^Bl 

0.03834 

130705 

2 . 11818 

29.7028 

4.86759 

1.27537 

0.16418 

0.00737 

ggjtHtlcld 

bbq 

■5^B1 

jg^BB 

bbi 

■iWiMilclH 

BBSS 

mximntw  i 

4.8673 

mm*Lw\ 

■lIKCEld 

-0.221 

0.0043 

-0.0512 

4.86692 

1.27385 

-0.2215 

■innwail 

-0.0601 

0 

-10 

25 

i 

1 

4.86684 

1.27366 

-0.0663 

0 

-20 

i 

1 

4.76611 

1.24259 

-0.1585 

-0 . 0097 

-0.0112 

0 

-20 

-12 

i 

1 

2.1249 

4.85204 

1.24019 

-0.1496 

0 

-20 

i 

1 

-0.0374 

0 

-20 

7 

i 

1 

-0.1411 

-0.0513 

0 

-20 

14 

i 

1 

130939 

2 . 12579 

-0.1483 

-0.0662 

0 

-20 

25 

i 

1 

131231 

2.13528 

1.2353 

0.0145 

-0.0785 

0 

-30 

i 

1 

131174 

2 . 13342 

4.7498 

1.2057 

0 

-30 

WBf^EEk 

IB^BI 

131140 

2.13232 

0.21414 

-0.0342 

-0.0108 

-0.009 

0 

-30 

H^BI 

^^bbi 

131182 

2.13368 

■iIWIcmJ 

-0.0376 

0 

-30 

3E^B1 

131136 

2.13219 

4.82815 

1.18613 

-0.0616 

0 

-30 

0.03844 

131068 

2 . 12997 

29.7853 

4.83239 

1.19572 

0.21345 

0.01023 

-0.0545 

-0.0846 

0 

-30 

25 

i 

i 

0.03843 

131033 

2 . 12884 

29.7774 

4.83563 

1.20305 

0.2143 

0.01219 

-0.0814 

0.0253 

-0.1057 

0 

14 

-20 

i 

i 

0.03849 

131219 

2.13489 

29.8197 

4.84565 

1.4226 

0.17478 

0.01054 

-0.4975 

0 . 01242 

-0.0761 

0 

14 

-12 

i 

i 

0.03855 

131414 

2 . 14123 

29.864 

4.9319 

1.42089 

0.17347 

0.00991 

-0.4969 

0 . 0087 

-0.0659 

0 

14 

0 

i 

i 

0.03852 

131331 

2.13851 

29.845 

4.93476 

1.42736 

0 . 1724 

0.00882 

-0.5031 

0.00093 

-0.0447 

0 

14 

7 

i 

i 

0.03847 

131143 

2.1324 

29.8023 

4.93432 

1.42638 

0.17154 

0.00799 

-0.4987 

-0 . 0035 

-0.0322 

0 

14 

14 

i 

i 

0.03842 

130997 

2 . 12766 

29.7692 

4.93298 

1.42333 

0.17013 

0.00749 

-0.4902 

-0.0068 

-0.0233 

0 

14 

25 

i 

i 

0.03845 

131079 

2 . 13031 

29.7877 

4.92834 

1.41284 

0.16789 

0.00705 

-0.4742 

-0 . 0096 

-0.0157 
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Table  19.  Matrix  Run  Data,  Tail  3 


Yaw 

5e 

5m 

Tail 

stabalizers 

M# 

Re# 

q_c 

Uoo 

alpha_c 

CL 

Cd 

C_, 

Cm 

Cn 

C_Y\r 

0 

0 

-30 

2 

0 

0.03939 

134835 

2.23741 

30.4136 

4.76786 

1.24656 

0.14607 

0.00591 

-0.34 

■mmr.-ictl 
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Table  20.  Matrix  Run  Data,  Tail  3  w/stabilizer 
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Appendix  B:  Balance  Calibration  Plots 


The  following  method  was  used  by  Mr.  Gehring  to  calibrate  and  install  the 
balance.  For  the  starting  point  of  this  process  it  is  assumed  that  the  balance  is  not 
attached  to  anything  and  the  sting  (balance  holder)  is  removed  from  the  tunnel.  These 
steps  are  not  meant  to  replace  proper  training  as  important  information  may  have  gone 
unobserved. 

Excitation  voltage  is  first  changed  within  the  tunnel  controller  program. 
Excitation  voltages  are  specific  for  the  balance  in  question,  with  the  AFIT-1  balance 
using  5V.  By  first  changing  the  voltage,  problems  with  forgetting  to  change  the  voltage 
later  are  nullified. 

Wires  attached  to  the  balance  must  first  be  threaded  through  the  balance 
holder.  These  wires  are  very  fine  and  to  not  have  enough  stiffness  on  their  own  to  push 
through  the  holder.  To  overcome  this  problem,  shrink  tubing,  roughly  6  inches  long,  was 
placed  over  the  ends  of  the  wire.  Care  had  to  be  taken  before  this  step  to  assure  the  wire 
marker  tags  were  slid  to  different  positions  so  they  would  not  bunch  and  fonn  a  blockage. 
After  using  the  shrink-tubed  wire  to  push  the  initial  length  of  wire  through  the  sting,  the 
shrink  tubing  was  removed  using  a  fine-pointed  pair  of  tweezers.  Tweezers  were  slid 
into  the  opening  of  the  shrink  tubing  and  pulled  upward  to  cut  it.  The  fine  point  would 
miss  any  wires  so  they  would  not  be  cut.  Upon  pulling  the  rest  of  the  balance  wire 
through  the  holder  the  balance  is  ready  to  be  attached  to  the  sting. 

Within  the  balance  are  a  number  of  screws  that  attach  the  balance  to  the  sting. 
These  screws  can  be  grouped  into  three  groups:  “push  in”,  “vertical”,  and  “push  out”. 
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The  names  refer  to  the  way  the  screws  and  their  holes  are  canted  and  how  they  appear  to 
affect  the  balance. 

While  inserting  the  balance  into  the  sting,  the  screw  holes  must  be  aligned 
properly.  The  alignment  of  these  holes  can  be  seen  in  the  right  edge  of  Figure  19.  When 
this  is  accomplished,  the  “push  in”  screws  are  finger  tightened  first.  Following  are  the 
vertical  screws  and  then  the  “push  out”  screws.  This  entire  assembly  is  then  slid  into  the 
calibration  holder,  located  adjacent  to  the  wind  tunnel  in  room  LI 54.  After  the  balance 
assembly  is  securely  fastened  into  the  calibration  fixture,  the  balance  wires  can  then  be 
attached  to  the  power  terminals.  Upon  applying  power  it  is  advisable  to  allow  24  hours 
for  the  internal  components  of  the  balance  to  warm  up. 

After  the  24  hour  period,  the  tunnel  controller  program  is  run  to  zero  out  the 
balance.  This  is  accomplished  by  changing  preset  voltages  in  the  block  diagram  portion 
of  the  program.  The  controller/data  acquisition  program  is  started  and  a  null  of  the 
balance  is  taken.  The  maximum  load  for  a  particular  direction  is  applied.  Voltage  gains, 
located  on  the  front  panel  of  the  controller/acquisition  program  are  changed  so  the  read¬ 
out  displays  the  proper  load.  The  load  is  then  removed  and  the  program  is  again  checked 
to  ensure  the  read-out  returns  to  zero.  This  is  repeated  for  each  of  the  other  force  and 
moment  directions  taken  by  the  balance. 

Calibration  of  the  balance  is  a  lengthy  process.  Pitch  itself  has  4  configurations  to 
be  tested.  The  first  is  orienting  the  balance  with  nonnal  positive  facing  up.  When  the 
pitch  calibration  fixture  is  applied  two  pitch  moments  are  calibrated,  positive  and 
negative.  The  balance  is  then  inverted  so  positive  nonnal  force  is  down.  The  pitch 
moment  is  then  calibrated  again.  Yaw  moment  also  has  4  configurations,  very  similar  to 
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pitch.  The  normal  force,  side  force,  and  roll  all  have  2  configurations.  Normal  direction 
is  with  the  positive  markings  on  the  balance  up  and  down.  Side  force  is  with  the  either 
side  of  the  balance  up  and  down.  Roll  is  calibrated  clockwise  and  counter  clockwise. 
Axial  force  is  tested  by  attaching  a  pulley  mechanism  to  the  balance  calibration  fixture. 

A  string  is  then  attached  to  the  balance  and  then  over  the  pulley.  Weights  are  then  added 
to  the  end  of  the  string  resulting  in  an  axial  force  applied  to  the  balance.  Graphs  were 
taken  of  the  calibration  before  the  MAV  tests.  The  graphs  in  Figures  69-79  demonstrate 
the  varying  directions  the  balance  is  calibrated. 


Positive  Normal  Force  (Lift) 


Figure  69.  Normal  Force  Calibration  Curve  (positive) 
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Negative  Normal  Force  (downforce) 


Figure  70.  Normal  Force  Calibration  Curve  (negative) 

Figure  69  and  Figure  70  were  produced  in  the  manner  described  previously.  The 
positive  normal  force  sensor  on  the  balance  was  positioned  facing  up  and  weights  hung 
from  it.  To  gather  the  negative  data,  the  positive  was  placed  down  and  weights  hung 
from  the  sensor.  The  remaining  graphs  from  this  calibration  are  included  below. 

These  plots  show  that  the  balance  and  software  gains  have  been  adjusted  properly 
to  give  an  indication  of  the  actual  load  applied  to  the  balance.  In  cases  where  a  spike  is 
present  in  the  data,  as  in  Figure  74  and  Figure  75,  this  is  from  a  number  of  smaller 
weights  being  removed  and  one  larger  weight  being  added.  Smaller  spikes  at  the 
beginning  of  each  step  in  weights  are  expected  as  the  weight  is  placed  on  the  weight  tray 
and  then  transient  responses  from  the  weight  addition  diminish. 
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Pounds  (5#'s  Max]  Pounds  (5#'s  max] 


Negative  Axial  Force  (Thrust) 


Figure  71.  Axial  force  Calibration  Plot 


Positive  Side  Force  (Right  Wing  Pull) 


Figure  72.  Positive  Side  Force  Calibration  Plot 
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inch-lbs  (10  in-ffs  max’  Pounds  (5#'s  Max] 


Negative  Side  Force  (Right  Wing  Push) 


Figure  73.  Negative  Side  Force  Calibration  Plot 


Positive  Pitch  Moment  (pitch  up) 


Figure  74.  Positive  Pitch  Moment  Calibration  Plot 
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Negative  Pitch  Moment  (pitch  down) 
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Figure  75.  Negative  Pitch  Moment  Calibration  Plot 


Positive  Yaw  Moment  (Nose  right) 


Figure  76.  Positive  Yaw  Moment  Calibration  Plot 
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inch-pounds  (4  in-#'s  Max)  inch-lbs  (5  in-#'s  Max 


Negative  Yaw  Moment  (Nose  Left) 


Figure  77.  Negative  Yaw  Moment  Calibration  Plot 


Positive  Roll  Moment  (Right  wing  down) 


Figure  78.  Positive  Roll  Moment  Calibration  Plot 
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inch-lbs  (4  in-#'s  Max) 


Negative  Roll  Moment  (right  wing  up) 


Figure  79.  Negative  Roll  Moment  Calibration  Plot 
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Appendix  C  Additional  Plots 


Plots  in  this  appendix  are  not  referenced  in  the  text.  These  plots  are  included  to 
better  show  stabilizer  effect  on  the  MAV/rotatable  tail  combination.  Figures  80  to  83 
show  similar  trends  as  discussed  for  Cl  plots  in  the  text.  Figure  83  departs  slightly  from 
trends  by  showing  an  increase  in  lift  at  5e  =  55°.  Change  in  tail  deflection  due  to  wind  is 
most  likely  the  cause. 


CL  due  to  Elevator  Deflection  (5e); 
20  mph  Tail  1  w /  stabilizers 


Negative  Elevator  Deflection,  5e,  (degrees) 
Figure  80.  CL  vs.  5C,  20  mph,  Tail  1  w/stabilizers 
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CL  due  to  Elevator  Deflection  (5e); 
20  mph  Tail  3  w /  stabilizers 


Figure  81.  CL  vs.  5e,  20  mph,  Tail  3  w/stabilizers 

CL  due  to  Elevator  Deflection  (5e); 

30  mph  Tail  1  w /  stabilizers 


Figure  82.  CL  vs.  5e,  30  mph,  Tail  1  w/stabilizers 
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CL  due  to  Elevator  Deflection  (5e); 
30  mph  Tail  3  w /  stabilizers 


Figure  84  through  87  are  the  Cd  plots  for  various  flight  conditions  and  tail 
configurations.  Drag  generally  increased  for  all  plots  at  5e  became  more  negative.  Figure 
86  and  Figure  87  again  show  a  decrease  in  drag  at  6e  =  55°.  This  corresponds  to  the 
increase  in  lift  seen  in  the  Cl  plots  for  the  same  condition. 
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CD  due  to  Elevator  Deflection  (5e); 
20  mph  Tail  1  w /  stabilizers 


Figure  84.  CD  vs.  5e,  20  mph,  Tail  1  w/stabilizers 


CD  due  to  Elevator  Deflection  (5e); 


Figure  85.  CD  vs.  5e,  30  mph,  Tail  1  w/stabilizers 
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CD  due  to  Elevator  Deflection  (5e); 
20  mph  Tail  3  w /  stabilizers 


Negative  Elevator  Deflection,  5e,  (degrees) 
Figure  86.  CD  vs.  5e,  20  mph,  Tail  3  w/stabilizers 


CD  due  to  Elevator  Deflection  (5e); 
30  mph  Tail  3  w /  stabilizers 


Negative  Elevator  Deflection,  5e,  (degrees) 
Figure  87.  CD  vs.  5e,  30  mph,  Tail  3  w/stabilizers 
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Figure  88  through  92  are  Cm  plots  for  the  various  configurations.  Figure  88 
shows  a  direct  comparison  of  one  flight  condition  for  both  configurations.  The  stabilizer 
weight  is  seen  to  reduce  the  pitching  moment  for  most  of  the  range  of  5e,  though  the  shift 
in  CG  was  somewhat  compensated  for  in  the  data  reduction  program.  A  sharp  drop  in  Cm 
is  seen  in  Figure  90  and  Figure  91  at  5e  =  55°.  This  anomaly  corresponds  with  similar 
anomalies  in  the  Cl  and  Cd  plots  above.  Similarly  the  anomaly  corresponds  with  trends 
in  the  text  that  were  related  to  a  flat  plate  at  incidence.  For  the  conditions  in  Figure  91 
the  largest  Cm  value  is  at  5e  =  45°  corresponding  with  the  largest  Cn0miai  seen  in  the  flat 
plate.  The  drop  off  after  5e  =  45°  in  Figure  91  is  similar  to  the  drop  in  Cn0nnai  at  the  same 
angular  condition. 


Cm  due  to  Elevator  Deflection  (5e); 
20  mph  Tail  1  w /  and  w/o  stabilizer 


Figure  88.  Cm  vs.  5e,  20  mph,  Tail  lw/&  w/o  stabilizer 
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Cm  due  to  Elevator  Deflection  (5e); 
30  mph  Tail  1  w /  stabilizers 


Negative  Elevator  Deflection,  5e,  (degrees) 
Figure  89.  Cm  vs.  5e,  30  mph,  Tail  1  w/stabilizers 


Cm  due  to  Elevator  Deflection  (5e); 
30  mph  Tail  1 


Negative  Elevator  Deflection,  5e,  (degrees) 
Figure  90.  Cm  vs.  5e,  30  mph,  Tail  1 
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Cm  due  to  Elevator  Deflection  (5e); 
20  mph  Tail  3  w/  stabilizers 


Figure  91.  Cm  vs.  5e,  20  mph,  Tail  3  w/stabilizers 


Cm  due  to  Elevator  Deflection  (5e); 
30  mph  Tail  3  w /  stabilizers 


Figure  92.  Cm  vs.  5e,  30  mph,  Tail  3  w/stabilizers 
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Directional  Stability  Comparison  of  Tail  1  vs  Tail  3 


Figure  93  shows  that  the  stabilizer  improved  stability.  Cnv|/  should  be  negative  for 
stability.  The  slopes  of  the  stabilizer  plots  in  this  figure  are  negative  while  without 
stabilizer  the  slope  is  positive. 
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Appendix  D  MATLAB®  Code  for  Data  Reduction 


0/^5i' 5jsjjc>jsjjs5jssic>}s5is5}s  |  -j-  ^  Cspt  Anthony  DgL/UCcI 

%*  Adapted  for  the  Balance  AFIT  1  and  rotary  tails  by  Lt.  Rivera  Parga  %*********** 
%*******  re-adapted  for  angled  balance  mount  by  ENS  Troy  Leveron  ********** 

0/^  Qp  J  j  pj-  Ivlomcnts 

0/^ h* h® h* h* h* h* h* h* h*  pj  WING  Prop  OFF'  AT ^PH^A.  S^VEEPS  H**i*H*H*H*5i*H*®i*H*5j*H*H*H*HsH*H*H**i* 


%This  Code  will  transfer  measured  Forces  and  Moments  on  the  AFIT  1  balance  to  Wind 

%(earth)  centered  frame  of  reference  by  correcting  for  tare  effects,  balance 

%interactions,  and  wind  tunnel  irregularities,  then  gives  a  file  with  all  the 

%corrected  data 

clear; 

clc; 

close  all; 
format  long 


%##################################################################### 

%  INPUT  DECK 

%FIRST  FILL  THE  FOLLOWING  INFORMATION 


Masskg=0.415;  %  Mass  of  the  UAV  in  KGS 

Troom  =  mean([72.3])  +  459.67;  %deg  R  **Changed  for  each  day  of  testing* 

P_barro  =  mean([29.06])  *  0.491 1541;  %Psi  ****Changed  for  each  day  of  testing**** 

%  Added  24  April  to  allow  easier  change  of  input  tare  and  data  fdes 

load  tarefile.txt;  %tarefile 

TareFile  =  tarefile(:,l:9); 

load  datafile.txt;  %datafile 

DataFile  =  datafile(:,l:9); 


%Offset  distances  from  the  Mounting  Block  to  the  Model  C.G.  (inches) 

Y_cmb  =  -(0.05);  %inches 
X_cmb  =  0.71;  %inches 
Z_cmb  =  -1.31;  %inches 

%  Requeried  for  the  Solid  body  blockage  corrections  due  to  wing  and  fuselage 
BodyVolume  =  ((9.42962435  *(2/1 6))+(0)+  (0.9375)+... 
(3.5*((2.75+1.35)/2)*((1.95+1.32)/2))... 

+(5. 25*3. 0*((1. 95+1. 45)/2)))  /  12A3; 

%(ftA3):  Tail+vertical  stabilizers  (Tail  2)+  Connector+Prop-to-Wing+Wing  Front-to- 
Wing  Back 

%  Requeried  for  the  Pitching  Moment  Correction 

I  t  =  9/12;  %  ft  =  length  from  tail  MAC  to  aircraft  CG 

Span  t  =(4+(6/16))  /  12;  %  ft  =  horizontal  span 

Tail_Area  =  (9.42962435)  /  144;  %  ftA2  =  horizontal  tail  area 

%  BEFORE  CONTINUING  IT  IS  NECCESARY  TO  CHANGE  THE  NAME: 

%  INPUT  DATA  FILE  AND  INPUT  DATA  TARE  FILE 
%  THE  OUTPUT  DATA  FILE 
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%###################################################################### 

%I.  Polynomial  Curve  fit  of  Chord  Dimensions  to  Spanwise  length  to  get  an 
%  Equation  defining  thet  Chord  anywhere  along  the  span. 

Span  Location  =  [0  2.62  (2.62*2)  (2.62*3)  (2.62*4)];  %in 
ChordDist  =  [1/2  59/16  5  89/16  6];  %in 
cl  =  polyfit(Span_Location,Chord_Dist,2); 

Chord  1  =  polyval(cl,Span_Location); 
c2  =  polyfit(Span_Location,Chord_Dist,3); 

Chord2  =  polyval(c2,Span_Location); 
c3  =  polyfit(Span_Location,Chord_Dist,4); 

Chord3  =  polyval(c3,Span_Location); 
figure(l) 

plot(Span_Location,Chord_Dist,'x',Span_Location, Chord  1,'-  .’,Span_Location,Chord2,'o- 
.',Span_Location,Chord3, ’*-'); 

%4th  order  chord  equation  as  a  function  of  the  span  (b)  or  C(y)  integrated 
%from  0  to  b/2  to  calculate  1/2  of  the  wing  area. 

chord_eqn  =  inline(’-0. 0004421 336783  l*b.A4+0. 01 737585355753  *b.A3- 
0.25 190354097469*b. A2+l. 765267 17557252*b+0. 50’); 

Wing_Area  =  (2  *  quad(chord_eqn, 0,4*2. 62))  /  144;  %ftA2 

%II.  Room  Conditions  and  Model  Specifics  : 

%  UNITS  are  in  Ft,  Sec,  lbm,  Psf,  Rankine,  fps 

Mass  =  (Masskg  *  1000)  *  0.0022046;  %lbm  (flex  MAV  with  batteries) 

Gas_Const=  1716;  %ft-lbf/Slug-R 

Density  =  (P_barro  *  144)/(1716  *  T_room);  %lbm/ftA3  or  lbf-sA2/ftA4 

Root  Chord  =  6  *  (1/12);  %ft 

Span  =  24  /  12;  %ft 

Aspect_Ratio  =  SpanA2  /  Wing_Area; 

Viscosity  =  .372e-6;  %slug/ft-s 

Speed  of  Sound  =  sqrt(1.4  *  Troom  *  Gas  Const);  %fps 

%III.  Solid  body  blockage  corrections  due  to  wing  and  fuselage 

K_1  =  0.9; 

K_3  =  0.93; 
delta  =  0.1 125; 

Taul  =0.83125; 

X_Section  =  (31/12)*(44/12);  %ftA2 

Wing  Volume  =  Wing_Area  *  (.006/12);  %ftA3 

Epsilon_sb_w  =  (K_l*Tau_l*Wing_Volume)  /  X_SectionA(3/2); 

Epsilon_sb_b  =  (K_3*Tau_l*Body_Volume)  /  X_SectionA(3/2); 

Epsilontot  =  Epsilonsbw  +  Epsilonsbb; 

%  III.  Load  the  static  tare  data  for  the  alpha  sweep  w/o  the  wind, 
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%  separate  each  force  from  the  file,  and  fit  a  4th  order  poly 
%  as  an  x-y  plot  (AoA  vs. Force)  for  each  of  the  6  force  sensors. 


%load  tarel  .txt;  %  Raw  tare  data  file  to  be  read  in;  commented  to  reflect 

change  above 

FILE=TareFile(:,l:9);  % 

j=i; 

k=l; 

L=length(FILE); 


for  i=l  :L  %Run  for  all  data  points  #  of  rows 

if  i~=L  %if  current  row  is  not  last  row,  go  to  next 

NEXT=i+l ;  %set  next  equal  to  the  value  of  the  next  row 

VALUE2=FILE(NEXT,  1 );  %set  value2  as  next  row  column  1 

else  if  i==L  %unless  the  it  is  the  last  value 

VALUE2=50;  %value2  set  to  50  to  end  the  sequence 

end 
end 

A(j,:)=FILE(i,:);  %set  row  j  of  A  equal  to  row  i  of  FILE 

VALUE  l=FILE(i,l);  %set  value  1  equal  to  row  i  column  1  of  FILE 

if  VALUE  1==VALUE2  %if  value  1  equals  value2,  go  to  next  row 

j=j+l; 

else  if  VALUE  1~=VALUE2  %if  value  1  and  value2  are  different  check 

if  length(A(:,  1))<5  %if  less  than  20  values,  ignored  due  to  angle  change 

j=i; 

clear  A; 

else  if  length(A(:,  1))>5  %if  more  than  20  values 

C=length(A(:,  1));  %find  length  of  A 

for  m=l  :9  %  Average  all  rows  of  the  like  values  in  A 

B(k,m)=mean(A(4:C,m));  %disregarding  first  10  for  vibrations 
end 

j=i; 

k=k+l; 
clear  A 
end 
end 

end 

end 

end 

if  B(k- 1 , 1  )<B((k-2),  1 ) 

B=B(l:(k-2),:) 

end 
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tare=[B]; 


% _ End  of  inserted  code 

[row, col]  =  size(tare);for  k  =  1  :row; 
theta_tare(k,:,:)  =  tare(k,  1  ).*  (pi/180); 

NF_tare(k,:,:)  =  tare(k,4); 

PM_tare(k,:,:)  =  tare(k,5); 

SF_tare(k,:,:)  =  tare(k,7); 

YM_tare(k,:,:)  =  tare(k,8); 

AF_tare(k,:,:)  =  tare(k,6); 

RM_tare(k,:,:)  =  tare(k,9); 
end 

NF_poly  =  polyfit(theta_tare,NF_tare,4); 

PM_poly  =  polyfit(theta_tare,PM_tare,4); 

SF_poly  =  polyfit(theta_tare,SF_tare,4); 

YM_poly  =  polyfit(theta_tare,YM_tare,4); 

AF_poly  =  polyfit(theta_tare,AF_tare,4); 

RMpoly  =  polyfit(theta_tare,RM_tare,4); 

%###################################################################### 

%IV.  Load  the  specific  test  run  files, 

clear  ('AA'/B'/C'/L') 

% _ 

%load  datal.txt; 

FILE=DataFile(:,:); 

j=i; 

k=1; 

L=length(FILE); 


%  Raw  data  file  to  be  read  in  commented 
%  to  reflect  change  made  above 

% 


for  i=1  :L 
if  i~=L 
NEXT=i+1 ; 


%Run  for  all  data  points  #  of  rows 
%if  current  row  is  not  last  row,  go  to  next 


%set  next  equal  to  the  value  of  the  next  row 
VALUE2=FILE(NEXT,1);  %set  value2  as  next  row  column  1;  changed  1  to  2  for  beta  runs 
else  if  i==L  %unless  the  it  is  the  last  value 

VALUE2=50; 


%value2  set  to  50  to  end  the  sequence 


end 
end 

A(j,:)=FILE(i,:); 

VALUE1  =FILE(i,1 ); 
if  VALUE  1—VALUE2 
j=j+1; 

else  if  VALUE  1~=VALUE2 
if  length(A(:,1))<5 


%set  row  j  of  A  equal  to  row  i  of  FILE 

%set  valuel  equal  to  row  i  column  1  of  FILE;  change  1  to  2  for  beta 
%if  valuel  equals  value2,  go  to  next  row 

%if  valuel  and  value2  are  different  check 
%if  less  than  20  values,  ignored  due  to  angle  change; 
%change  to  2  for  beta 


j=i; 

clear  A; 

else  if  length(A(:,1))>5  %if  more  than  20  values;  change  1  to  2  for  beta  runs 

C=length(A(:,1));  %find  length  of  A;  change  1  to  2  for  beta  runs 

for  m=1 :9  %Average  all  rows  of  the  like  values  in  A 

B(k,m)=mean(A(4:C,m));  %disregarding  first  10  for  vibrations 
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end 

j=i; 

k=k+1 ; 
clear  A 
end 
end 
end 
end 
end 

if  B(k-1 ,1)<B((k-2),1)  %change  1  in  ",1)u  to  2  for  beta  runs 
B=B(1  :(k-2),:) 
end 

sample_data=[B]; 


% _ End  of  inserted  code 

% _ End  of  inserted  code 

[row2,col2]  =  size(sample_data); 
for  i  =  1  :row2; 

%Angles  of  the  model  during  test  runs  (Roll,  Pitch  {AoA},  Yaw  {Beta}): 
phi  =  0; 

theta(i,:)  =  sample_data(i,l)  .*  (pi/180);  %radians 

si(i,:)  =  sample_data(i,2)  .*  (pi/180);  %radians 

Wind_Speed(i,:)  =  sample_data(i,3)  .*  (5280/3600);  %fps 

%Flight  Parameters  (Re#,  Ma#,  Dynamic  Pressure): 

q  =  (.5  *  Density)  .*  Wind_Speed.A2;  %lbf/ftA2 

q_Corrected  =  q  .*  (1  +  Epsilon_tot)A2;  %lbf/ftA2 

Wind_Speed_Corrected  =  Wind_Speed  .*  (1  +  Epsilon_tot);  %fps 

Mach_Number  =  Wind_Speed_Corrected  ./  Speed_of_Sound;  %NonDimensional 

Reynolds  Number  =  ((Density  *  Root  Chord)  .*  Wind  Speed  Corrected)  ./Viscosity; 

%N  onD  imensional 

Flight_Parameters  =  [Mach_Number  Reynolds_Number  q_Corrected]; 

%individual  forces  and  moments  for  each  sensor: 

NF_test(i,:,:)  =  sample_data(i,4); 

PM_test(i,:,:)  =  sample_data(i,5); 

SF_test(i,:,:)  =  sample_data(i,7); 

YM_test(i,:,:)  =  sample_data(i,8); 

AF_test(i,:,:)  =  sample_data(i,6); 

RM_test(i,:,:)  =  sample_data(i,9); 

%###################################################################### 

%V.  Subtract  the  effect  of  the  static 
%  weight  with  the  tare  polynomials  above 

%Evaluating  the  actual  test  theta  angle  (AoA)  in  the  tare  polynomial  to 
%determine  the  tare  values  for  the  angles  tested  in  each  run. 

NF_eval  =  polyval(NF_poly, theta); 

PM_eval  =  polyval(PM_poly, theta); 
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SF_eval  =  polyval(SF_poly, theta); 

YM_eval  =  polyval(YM_poly,  theta); 

AF_eval  =  polyval(AF_poly,  theta); 

RM_eval  =  polyval(RM_poly,  theta); 

%The  Time-Averaged  (raw)  forces  and  momentums  NF,AF,SF,PM,YM  AND  RM 
measured  in  the  wind 

%tunnel  (body  axis)  with  the  tare  effect  of  the  weight  subtracted  off. 

NF_resolved  =  NF_test  -  (NF_eval); 

PM_resolved  =  PM_test  -  (PM_eval); 

SF_resolved  =  SF_test  -  (SF_eval); 

YM_resolved  =  YM_test  -  (YM_eval); 

AF_resolved  =  AF_test  -  (AF_eval); 

RM_resolved  =  RM_test  -  (RM_eval); 

Forces_minus_tare  =  [NF_resolved,  AF_resolved,  PM_resolved,  RM_resolved, 
YM_resolved,  SF_resolved]'; 

%VI.  CORRECT  FORCES  AND  MOMENTS  FOR  BALANCE  INTERATIONS  (body 
axis) 

%USING  THE  REDUCTION  EQUATIONS 

%LET  US  SET  A  MAXIMUN  NUMBER  OF  INTERATIONS  (FOR  AVOIDING  AN 
INFINIT  LOOP) 

MAXIT=100; 

%SET  THE  LIMIT  FOR  THE  DIFFERENCE  BETWEEN  INTERATIONS(CRITERIA 
FOR  FINISH  THE  INTERATIONS) 

LIMIT=  10E-14; 

%MATCHING  EACH  NAME  WITH  THE  DATA 
MNF=NF_resolved(i); 

MAF=AF_resolved(i); 

MPM=PM_resolved(i); 

MRM=RM_resolved(i) ; 

MYM=YM_resolved(i); 

MSF=SF_resolved(i); 

%INPUT  OF  THE  CONSTANTS  VALUES  FROM  THE  MATRIX  FOR 
SENSITIVITIES  AND  INTERATIONS 

K=[0  -1.3567E-03  -3.8021E-03  -4.2814E-03  -1.6966E-03  1.7567E-03  ... 

5.3167E-05  -1.3867E-04  -5.5629E-05  3.5181E-05  1.0601E-05  -2.5271E-04... 
5.6693E-05  -1.9537E-04  1.7908E-05  -3.6606E-05  -4.9934E-05  4.1205E-05... 
2.5648E-05  -1.9289E-05  8.9661E-05  -1.9594E-05  -4.9859E-04  -1.1599E-03... 
5.7163E-05  8.9798E-05  -7.8591E-05  9.3187E-03  0  -3.8421E-03  3.5740E-03... 
9.7714E-05  -2.7776E-03  -1.3552E-04  5.1538E-04  2.2082E-04  -1.2706E-05... 
-2.3637E-05  1.3686E-05  1.1085E-04  -3.6557E-06  4.9876E-06  8.1085E-06... 

3.7381E-05  1.2791E-04  -9.4527E-06  -2.3083E-06  -1.2046E-06  7.8161E-04... 
-1.1997E-03  -3.0560E-05  -6.6202E-05  3.7227E-04  -2.1469E-04  4.8386E-03... 
-3.7387E-03  0  -1.8479E-02  3.9077E-03  9.9165E-04  -1.4825E-05  -1.4830E-06... 
6.0845E-05  8.0667E-05  1.8547E-05  -5.0212E-05  1.0539E-04  -2.2676E-04... 
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4.3793E-05  -1.0456E-05  -8.1186E-06  -2.1653E-05  -3.3070E-05  1.7280E-05... 
-7.4509E-05  -3.4399E-05  -8.2999E-04  -6.7962E-04  4.052 IE-05  -5.1604E-05... 

9.1 132E-06  -5.7360E-03  -2.2213E-04  9.9131E-04  0  -9.5790E-03  6.71 14E-03... 
3.6824E-05  1.0056E-04 -3.7105E-05  -9.0295E-05  -7.4580E-05  1.4814E-04... 

7.2634E-05  -8.4778E-06  6.3486E-05  5.6328E-05  -1.3617E-04  2.2196E-05... 

1.3606E-05  -3.6689E-05  8.3283E-05  1.1865E-04  1.8544E-05  -1.9831E-05... 

1.7894E-05  -6.8164E-05  -7.0892E-05  1.2378E-03  1.6961E-03  -6.5102E-03... 
-9.3202E-03  0  5.1349E-03  1.3612E-05  -1.3175E-04  7.2442E-06  5.6705E-04... 
-1.4723E-05  -4.8656E-05  -1.4282E-04  5.971  IE-05  5.9046E-05  -3.6490E-04... 

7.488 IE-05  5.4601E-06  1.0129E-03  -1.3867E-04  8.1617E-05  6.6053E-05... 

-1.3417E-05  9.0025E-05  -4.5362E-05  -4.4672E-06  9.5087E-05  -3.4077E-02... 
7.9142E-04  1.6667E-03  -6.6512E-03  8.1538E-03  0  -1.4185E-05  7.3209E-05... 
-2.5849E-05  1.2325E-03  -4.1696E-05  4.6266E-05  8.6146E-05  2.1436E-05... 

5.0874E-05  -3.2738E-04  2.2218E-04  8.6478E-06  7.3395E-04  -4.1453E-05... 

3.5719E-05  2.5313E-05  1.5182E-04  3.6007E-05  -2.8844E-05  8.9741E-05... 
-7.3257E-05]; 

%COMPUTE  THE  UNCORRECTED  FORCES  AND  MOMENTS  BY 
%CONSIDERING  THAT  THE  PRIME  SENSITIVITY  CONSTANTS  ARE  ALREADY 
APLIED: 

NF1=MNF; 

AF1=MAF; 

PM1=MPM; 

RM1=MRM; 

195 

YM1=MYM; 

SF1=MSF; 

%FOR  THE  FIRST  INTERACTION  LET  US  INIZIALICE  THE  VALUES  OF  FORCES 
AND 

%MOMENTS  WITH  THE  VALUES  OF  THE  UNCORRECTED  FORCES  AND 

MOMENTS 

NF(1)=NF1; 

AF(1)=AF1; 

PM(1)=PM1; 

RM(1)=RM1; 

YM(1)=YM1; 

SF(1)=SF1; 

%DOING  THE  INTERACTION  EQUATIONS: 
for  n=2:MAXIT; 

NF(n)=NF  1  -((K(2)*AF(n- 1  ))+(K(3)*PM(n- 1  ))+(K(4)*RM(n- 1  ))+(K(5)*YM(n- 
1  ))+(K(6)*SF(n- 1  ))+(K(7)*NF(n- 1  )A2)+. . . 

(K(8)*(NF(n-l)*AF(n-l)))+(K(9)*(NF(n-l)*PM(n-l)))+(K(10)!l=(NF(n-l)*RM(n- 
1  )))+(K(  11)*  (NF(n-  l)*YM(n-l )))+ . . . 

(K(  1 2)  *  (NF  (n- 1 )  *  SF(n- 1  )))+(K(  1 3)  *  ( AF(n- 1 )  A2))+(K(  1 4)*  (AF  (n- 1 )  *PM(n- 
1  )))+(K(  1 5)*(AF(n- 1  )*RM(n- 1 )))+.. . 

(K(  1 6)  *  (AF  (n- 1 )  *  YM(n- 1  )))+(K(  1 7)  *  (AF  (n- 1 )  *  SF(n- 1  )))+(K(  1 8)  *  (PM(n- 
1  )A2))+(K(  1 9)  *  (PM(n- 1 )  *RM(n- 1 )))+. . . 
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(K(20)*(PM(n- 1  )*  YM(n- 1  )))+(K(2 1  )*(PM(n- 1  )*SF(n-  l)))+(K(22)*(RM(n- 
1  )A2))+(K(23 )  *  (RM(n- 1  )*YM(n- 1 )))+... 

(K(24)*(RM(n- 1  )*SF(n- 1  )))+(K(25)*(YM(n- 1  )A2))+(K(26)*(YM(n- 1  )*SF(n- 
l)))+(K(27)*(SF(n-l)A2))); 

AF(n)=AF  1  -((K(28)*NF(n- 1  ))+(K(30)*PM(n- 1  ))+(K(3 1  )*RM(n- 1  ))+(K(32)*  YM(n- 
l))+(K(33)*SF(n-l))+(K(34)  *NF  (n- 1  )A2 )+ . . . 

(K(35)*(NF(n-l)*AF(n-l)))+(K(36)*(NF(n-l)*PM(n-l)))+(K(37)*(NF(n- 
1  )*RM(n- 1  )))+(K(3  8 )  *  (NF  (n- 1  )*YM(n- 1 )))+... 

(K(39)*(NF(n-l)*SF(n-l)))+(K(40)*(AF(n-l)A2))+(K(41)*(AF(n-l)*PM(n- 
1  )))+(K(42)  *  ( AF(n- 1 )  *  RM(n- 1 )))+... 

(K(43 )  *  (AF  (n- 1 )  *  YM(n- 1  )))+(K(44)  *  (AF(n- 1 )  *  SF(n- 1  )))+(K(45)  *(PM(n- 
1  )A2))+(K(46)  *  (PM(n- 1  )*RM(n- 1 )))+... 

(K(47)*(PM(n- 1  )* YM(n-  l)))+(K(48)*(PM(n- 1  )*SF(n- 1  )))+(K(49)*(RM(n- 
l)A2))+(K(50)*(RM(n-l)*YM(n-l )))+... 

(K(51)*(RM(n-l)*SF(n-l)))+(K(52)*(YM(n-l)A2))+(K(53)*(YM(n-l)*SF(n- 

l)))+(K(54)*(SF(n-l)A2))); 

PM(n)=PM  1  -((K(55)*NF(n- 1  ))+(K(56)*  AF(n-  l))+(K(58)*RM(n- 1  ))+(K(59)*YM(n- 
1  ))+(K(60)*  SF(n- 1  ))+(K(6 1 )  *NF(n- 1  )A2)+. . . 

(K(62)  *(NF  (n- 1 )  *  AF  (n- 1  )))+(K(63  )*(NF  (n- 1 )  *PM(n- 1  )))+(K(64)  *  (NF(n- 
1  )*RM(n- 1  )))+(K(65)*(NF(n- 1  )*YM(n-l )))+... 

(K(66)*(NF(n4)*SF(n-l)))+(K(67)*(AF(n-l)A2))+(K(68)*(AF(n-l)*PM(n- 
1  )))+(K(69)*(AF(n-l)*RM(n- 1 )))+... 

(K(70)*(AF(n-  l)*YM(n- 1  )))+(K(7 1  )*(AF(n- 1  )*SF(n- 1  )))+(K(72)*(PM(n- 
1  )A2))+(K(73)*(PM(n- 1  )*RM(n- 1 )))+. . . 

(K(74)*(PM(n- 1  )*' YM(n-  l)))+(K(75)*(PM(n- 1  )*SF(n-  l)))+(K(76)*(RM(n- 
1  )A2))+(K(77)*(RM(n- 1  )*YM(n- 1 )))+... 

(K(78)*(RM(n-l)*SF(n-l)))+(K(79)*(YM(n-l)A2))+(K(80)*(YM(n-l)*SF(n- 

l)))+(K(81)*(SF(n-l)A2))); 

RM(n)=RM  1  -((K(82)*NF(n-l))+(K(83)*AF(n- 1  ))+(K(84)*PM(n-l))+(K(86)*YM(n- 
1  ))+(K(  8  7)  *  SF(n- 1 )  )+(K(  8  8 )  *NF  (n- 1  )A2 )+ . . . 

(K(89)*(NF(n- 1  )*  AF(n-  l)))+(K(90)*(NF(n- 1  )*PM(n-  1)))+(K(9  l)*(NF(n- 
1  )*RM(n- 1  )))+(K(92 )  *  (NF  (n- 1  )*YM(n- 1 )))+... 

(K(93)*(NF(n-l)*SF(n-l)))+(K(94)*(AF(n-l)A2))+(K(95)*(AF(n-l)*PM(n- 
1  )))+(K(96)*(AF(n- 1  )*RM(n- 1 )))+.. . 

(K(97)*(AF(n-  l)*YM(n- 1  )))+(K(98)*(AF(n- 1  )*SF(n- 1  )))+(K(99)*(PM(n- 
1  )A2))+(K(  1 00)*(PM(n- 1  )*RM(n- 1 )))+. . . 

(K(  1 0 1  )*(PM(n-  l)*YM(n- 1  )))+(K(  1 02)*(PM(n- 1  )*  SF(n- 1  )))+(K(  1 03)*(RM(n- 
1  )A2))+(K(  1 04 )  *  (RM(n-  l)*YM(n- 1 )))+... 

(K(  1 05)  *  (RM(n- 1 )  *  SF(n- 1  )))+(K(  1 06)  *  ( YM(n- 1  )A2))+(K(  1 07)  *  (YM(n- 1 )  *  SF  (n- 
l)))+(K(108)*(SF(n-l)A2))); 

YM(n)= YM 1  -((K(  1 09)*NF  (n- 1  ))+(K(  1 1 0)  *  AF(n- 1  ))+(K(  1 1  l)*PM(n- 
1  ))+(K(  1 1 2)*RM(n- 1  ))+(K(  1 1 4)*  SF(n- 1  ))+(K(  1 1 5)*NF(n- 1  )A2)+. . . 

(K(l  16)*(NF(n-l)*AF(n-l)))+(K(l  17)*(NF(n-l)!l!PM(n-l)))+(K(l  18)*(NF(n- 
1  )*RM(n- 1  )))+(K(  1 1 9)*(NF(n-  l)*YM(n- 1 )))+.. . 

(K(120)*(NF(n-l)*SF(n-l)))+(K(121)*(AF(n-l)A2))+(K(122)*(AF(n-l)*PM(n- 
1  )))+(K(  123)*(AF(n-l)*RM(n-l )))+ . . . 
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(K(  1 24) *  *  ( AF  (n- 1 )  *  YM(n- 1  )))+(K(  125)*  ( AF  (n- 1 )  *  SF(n- 1  )))+(K(  1 2  6)  *  (PM(n- 
1  )A2))+(K(  127)*(PM(n- 1  )*RM(n- 1 )))+. . . 

(K(128)*(PM(n-l)*YM(n-l)))+(K(129)*(PM(n-l)*SF(n-l)))+(K(130)*(RM(n- 
1  )A2))+(K(  131)*  (RM(n- 1 )  *  YM(n- 1 )))+.. . 

(K(  1 3  2)  *  (RM(n-  1 )  *  SF(n- 1  )))+(K(  1 3  3  )  *  ( YM(n-  1)A2))+(K(134)*  ( YM(n- 1 ) *  SF(n- 
1  )))+(K(  1 3  5)  *  (SF(n- 1  )A2))); 

SF(n)=SF  1  -((K(  136)  *NF  (n- 1  ))+(K(  1 3  7)  *  AF  (n- 1  ))+(K(  1 3  8)  *PM(n- 1  ))+(K(  1 3  9)  *RM(n- 
1  ))+(K(  1 40)  *  YM(n- 1  ))+(K(  1 42)  *NF(n- 1  )A2)+ . . . 

(K(  143)*  (NF  (n- 1 )  *  AF  (n- 1  )))+(K(  1 44)  *  (NF  (n- 1 )  *  PM(n- 1  )))+(K(  1 45  )  *  (NF(n- 
1  )*RM(n- 1  )))+(K(  1 46)*(NF(n- 1  )*  YM(n- 1 )))+.. . 

(K(147)*(NF(n-l)*SF(n-l)))+(K(148)*(AF(n-l)A2))+(K(149)*(AF(n-l)*PM(n- 
1  )))+(K(  1 5  0)  *  ( AF  (n- 1 )  *  RM(n- 1 )))+.. . 

(K(151)*(AF(n-l)*YM(n-l)))+(K(152)*(AF(n-l)*SF(n-l)))+(K(153)*(PM(n- 
1  )A2))+(K(  1 54)*(PM(n- 1  )*RM(n- 1 )))+. . . 

(K(  1 55)*(PM(n- 1)*  YM(n- 1  )))+(K(  1 56)*(PM(n-  l)*SF(n-  1)))+(K(  1 57)*(RM(n- 
1  )A2))+(K(  1 5  8)  *  (RM(n- 1  )*YM(n- 1 )))+.. . 

(K(  1 5  9)  *  (RM(n- 1 )  *  SF(n- 1  )))+(K(  1 60)  *  ( YM(n- 1  )A2))+(K(  161)*  (YM(n- 1 )  *  SF  (n- 
1  )))+(K(  1 62)  *  (SF(n- 1  )A2))); 

%  SET  THE  LIMIT  FOR  THE  DIFFERENCE  BETWEEN  INTERATIONS(CRITERIA 
FOR  FINISH  THE  INTERATIONS) 

DIFFNF  (n)=abs(NF  (n)-NF  (n- 1 )) ; 

DIFFAF(n)=abs(AF(n)-AF(n- 1)); 

DIFFPM(n)=abs(PM(n)-PM(n- 1 )); 

DIFFRM(n)=abs(RM(n)-RM(n- 1)); 

DIFFYM(n)=abs(YM(n)-YM(n- 1 )); 

DIFFSF(n)=abs(SF(n)-SF(n-l)); 

if  DIFFNF(n)&DIFFAF(n)&DIFFPM(n)&DIFFRM(n)&DIFFYM(n)&DIFFSF(n)  < 

LIMIT 

break 

end 

end 

%disp('THE  FINAL  VALUES  ARE  (NF,AF,PM,RM,YM,SF):') 

Corrected_Data( :  ,i)=  [NF  (n) ;  AF  (n)  ;PM(n)  ;RM(n) ;  YM(n) ;  SF(n)  ] ; 

%disp('THE  FINAL  DIFFERENCE  BETWEEN  INTERATIONS  ARE(FOR 
%NF,AF,PM,RM,YM,SF)  :') 

%FINAL_DIFFERENCE=[DIFFNF(n),DIFFAF(n),DIFFPM(n),DIFFRM(n),DIFFYM(n), 

%DIFFSF(n)] 

%disp('THE  NUMBER  OF  INTERATIONS  USED  WAS:’) 

%VII.  Calculation  of  the  Axial,  Side,  &  Normal  Forces  from  the  corrected  balance 
%  forces  in  the  Body  Axis  reference  frame 

%###################################################################### 

Forces_b(:,i)  =  [Corrected_Data(2,i);  Corrected_Data(6,i);  Corrected_Data(l,i)]; 
%Calculation  of  the  Drag,  Side,  &  Lift  Forces  in  the  Wind  Axis  reference  frame 
Forces_w  =  [Forces_b(l,:).*cos(theta').*cos(si')+Forces_b(2,:).*sin(si’)+Forces_b(3,:). 

* sin( theta') . * cos(si');  -F orces_b(  1 , :).  * sin(si') . * cos(theta')+ 

Forces_b(2 , :) . * cos(si')  -Forces_b(3 , :) . * sin(theta') . * sin(si'); 
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-F  orces_b(  1 , *  sin(theta')+Forces_b(3 , .  *  cos(theta')] ; 

%First  entry  is  the  moments  calculated  by  the  balance  or  direct  calculation 
%in  the  Body  Reference  Frame.  Balance  measures  Roll  (1),  Yaw  is  about  the 
%z-axis  (n),  and  Pitch  is  about  the  y-axis  (m).  Distances  from  strain 
%gages  to  C.G.  are  in  INCFIES.  Moments  are  in-lbf 
m  =  Corrected_Data(3,i); 
n  =  Corrected_Data(5,i); 

1  =  Corrected_Data(4,i); 

Moments_b(:,i)  =  [1;  m;  n]; 

%Second  entry  is  the  conversion  from  the  "Balance  Centeric"  moments  to  the 
%Wind  Reference  monments  with  respect  to  the  Balance  Center  (be) 

Moments_w_bc  =  [Moments_b(l,:).*cos(theta').*cos(si’)-Moments_b(2,:). 
*sin(si')+Moments_b(3,:).*sin(  theta').  *cos(si'); 

Moments_b(  1 , :) .  *  sin(si') .  *  cos(theta')+Moments_b(2 , :) .  *  cos(si')+ 

Moments_b(3 , :) .  *  sin(theta').  *  sin(si'); 

-Moments_b(  1 , :) .  *  sin(theta')+Moments_b(3 , :) .  *  cos(theta')] ; 

%Finally,  the  balance  centered  moments  are  converted  to  moments  about  the 
%Model’s  Center  of  Mass  (cm)  or  Center  of  Gravity  (CG) 
cgdist=sqrt((X_cmb)A2+(Z_cmb)A2);  %Obtaining  the  direct  distance  between  the 
center  of  the  balance  and  the  center  of  mass 

w=atan(-Z_cmb/X_cmb);  %Obtaining  the  angle  between  cgdist  and  the  x  axes  at  zero 
angle  of  attack 

X_cm(i,:)=  cos(theta(i,:)+w)*cos(si(i,:))*(cgdist); 

Y_cm(i,:)  =  Y_cmb  +  X_cm(i,:)*tan(si(i,:)); 

Z_cm(i,  :)=  -sin(theta(i,  :)+w)  *  (cgdist); 

Moments_w_cg_u  =  [Moments_w_bc(l  ,:)  +  Z_  cm(i,:)*Forces_w(2,:)  +  Forces_w(3,:)* 
Y_cm; 

Moments_w_bc(2,:)  -  Forces_w(3,:)*  X_cm(i,:)  +  Forces_w(l,:)* 

Z_cm(i,:); 

Moments_w_bc(3,:)  -  Forces_w(l,:)*  Y_cm  -  Forces_w(2,:)* 

X_cm(i,:)]; 

%###################################################################### 

%VIII.  Calculation  of  the  actual  Lift  and  Drag  no  dimensional  Coefficients,  uncorrected 
for  tunnel  effects,  (Cl  and  Cd) 

C_D_u  =  Forces_w(l,:)  ./  (q_Corrected'  .*  Wing_Area); 

C_Y_u  =  Forces_w(2,:)  ./  (q_Corrected'  .*  Wing_Area); 

C_L_u  =  Forces_w(3,:)  ./  (q_Corrected'  .*  Wing_Area); 

Coefficients  =  [C_L_u;  C_D_u;  C_Y_u]’; 

Ave_Cl  =  mean(Coefficients(:,l)); 

Ave_Cd  =  mean(Coefficients(:,2)); 
end 

%IX  Drag  Coefficient  Correction 
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C_D_o  =  min(Coefficients(:,2)); 

C_L_u_sqrd  =  Coefficients(:,l).A2; 

Delta_C_D_w  =  ((delta  *  Wing_Area)  /  X_Section)  .*  C_L_u_sqrd; 

CDCorrected  =  C_D_u'  +  DeltaCDw; 

%X.  Angle  of  Attack  due  to  upwash  Correction 

alphabefore  =  sample_data(:,l); 

alpha  =  [alpha_before]-[5]  %changed  18  Apr  for  mounting  block  angle 
Delta_alpha_w  =  ((delta  *  Wing_Area)  /  X_Section)  .*  (57.3  *  C_L_u); 
alphaCorrected  =  alpha  +  Deltaalphaw'; 

%XI.  Pitching  Moment  Correction 

tau2  =  0.65; 

c_bar  =  (mean([6,  5+9/16,  5,  3+11/16,  0.5]))  /  12;  %  ft  =  Mean  Chord  of  wing 

Vbar  =  (Tail  Area  *  I  t)  /  (Wing  Area  *  c  bar);  %  Horizontal  tail  volume  ratio 

etat  =  1.0; 

epsilono  =  0; 

i_t  =  pi/4;  %  radians 

i_w  =  0; 

Aspect_Ratio_t  =  Span_tA2  /  Tail_Area; 

DepslionDalpha  =  ((2  .*  C_L_u)  ./  (pi*  AspectRatio))'; 
epsilon  =  epsilon  o  +  (D  epslion  D  alpha  .*  alpha  Corrected  ); 
alpha  t  =  alpha  Corrected  -  i_w  -  epsilon  +  i  t; 

C_L_alpha_t  =  ((0.1*  Aspect_Ratio)  /  (Aspect_Ratio_t  +2))  *  0.8; 
D_Cm_cg_t_D_alpha_t  =  -C  L  alpha  t*  V  bar  *  eta  t; 

Delta  C  m  cg  t  =  ((D_Cm_cg_t_D_alpha_t)  *  (delta*tau2)  *  (Wing  Area  / 
X_Section)  .*  (C_L_u  *  57.3))'; 

Cl_w_cg  =  Moments_w_cg_u(l,:)  ./  (q_Corrected’  .*  (Wing_Area  *  Span*  12)); 
Cm_w_cg_u  =  Moments_w_cg_u(2,:)  ./  (q_Corrected’  .*  (Wing_Area  *  c_bar*12)); 
Cn  w  cg  =  Moments_w_cg_u(3,:)  ./  (q_Corrected’  .*  (Wing_Area  *  Span*  12)); 

Cm  w  cg  corrected  =  Cm  w  cg  u  -  Delta  C  m  cg  t' ; 

Corrected  Moment  Coefficients  =  [Cl  w  cg'  Cm  w  cg  corrected'  Cn  w  cg’]; 
%OBTAINING  THE  MOMENTS  COEFFICIENTS  CORRECTED  ABOUT  THE 
CENTER  OF  THE 
%BALANCE 

Cl  w  bc  =  Moments_w_bc(l,:)  ./  (q_Corrected'  .*  (Wing_Area  *  Span*  12)); 
Cm_w_bc_u  =  Moments_w_bc(2,:)  ./  (q_Corrected'  .*  (Wing_Area  *  c_bar*12)); 
Cn  w  bc  =  Moments_w_bc(3,:)  ./  (q_Corrected'  .*  (Wing_Area  *  Span*  12)); 

Cm  w  bc  corrected  =  Cm  w  bc  u  -  Delta  C  m  cg  t' ; 

Corrected  Moment  Coefficients  bc  =  [Cl  w  bc'  Cm  w  bc  corrected’  Cn  w  bc’]; 
%XII.  OUTPUT  VARIABLES  FORMATING 
alpha  =  sample_data(:,1); 
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fprintf('  Mach  Number  Reynolds  Number  Dynamic  Pressure(Psf)\r') 

Flight_Parameters 
fprintf('  \r'); 

fprintf('  Loads  are  in  Ibf  and  arranged  [D  S  L]  across  the  top  and  increments  of  alpha  down  the 
side  \r') 

Forces_w' 
fprintf('  \r') 

fprintf('  Moments  are  in  in-lbf  and  arranged  [L  M  N]  down  the  side  and  increments  of  alpha  along 
the  top  \r') 

Moments_w_cg_u 
fprintfO  \r') 

fprintf('  Cl_u  Cd_u  CY_u  \r'); 

Coefficients 
fprintf('  \r') 

fprintf('  Del_CD_w  CD_u  CD_Corrected  \r'); 

Compare_CD  =  [Delta_C_D_w  C_D_u'  C_D_Corrected] 
fprintf('  \r') 

fprintf('  Del_alpha_w  alpha_g  alpha_Corrected  \r'); 

Compare_alpha  =  [Delta_alpha_w'  alpha  alpha_Corrected  ] 
fprintf('  \r') 

fprintf('  Cl_cg_wind  Cm_cg_corrected_w  Cn_cg_wind  \r'); 

Corrected_Moment_Coefficients 
fprintf('  \r') 

fprintf('  M#  Re#  q_c  Uoo  alpha_c  C_L  C_D_c  Cl_cg_w 
Cm_cg_c_w  Cn_cg_w  C_Y\r'); 

YY=[Flight_Parameters  (Wind_Speed_Corrected  .*  (3600/5280))  alpha_Corrected  C_L_u' 
C_D_Corrected  Corrected_Moment_Coefficients  C_Y_u']%pressure] 

%XX=['M#'  'Re#'  'q_c'  'Uoo'  'alpha_c'  'C_L'  'C_D_c'  'Cl_cg_w'  'Cm_cg_c_w'  'Cn_cg_w  \r']; 

%ZZ=[XX;  YY]; 

wkl  write('output.xls', YY, 2,0) 


Max_CI  =  max(Coefficients(:,1)); 

%  LET  US  SAVE  TOTAL  DATA  IN  A  EXTERNAL  FILE 
%dlmwrite('TEST  TOTAL  DATA  T2  DE  M9  DRN  20',TOTAL_DATA,'\t') 
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